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Abstract 

Conditions under which compression occurs and collective expansion devel- 
ops in energetic symmetric reactions of heavy nuclei, are analyzed, together 
with their effects on emitted light baryons and pions. Within transport sim- 
ulations, it is shown that shock fronts perpendicular to beam axis form in 
head-on reactions. The fronts separate hot compressed matter from normal 
and propagate into projectile and target. As impact parameter increases, 
the angle of inclination of the fronts relative to beam axis decreases, and 
in-between the fronts a weak tangential discontinuity develops. Hot matter 
exposed to the vacuum in directions perpendicular to shock motion (and par- 
allel to fronts), starts to expand sideways, early within reactions. Expansion 
in the direction of shock motion follows after the shocks propagate through 
nuclei, but due to the delay does not acquire same strength. Expansion af- 
fects angular distributions, mean-energy components, shapes of spectra and 



-mail: danielewicz@nscl.nscl.msu.edu 



1 



mean energies of different particles emitted into any one direction, and fur- 
ther particle yields. Both the expansion and a collective motion associated 
with the weak discontinuity, affect the magnitude of sideward flow within re- 
action plane. Differences in mean particle energy components in and out of 
the reaction plane in semicentral collisions, depend sensitively on the rela- 
tive magnitude of shock speed in normal matter and speed of sound in hot 
matter. The missing energy, considered in the past in association with low 
pion-multiplicity in central reactions, may be identified with the energy of 
collective expansion. Relations are established which approximately govern 
the behavior of density and entropy in the compressed region in reactions with 

beam energy and impact parameter. 
PACS numbers: 27.75.+r, 21.65.+f, 24.10.-i 
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I. INTRODUCTION 



Within the past two decades a considerable effort in the area of energetic heavy-ion re- 
actions was put into identification of signatures of collective phenomena in excited nuclear 
matter. Particular effort was devoted to the effects of collective motion. The latter has 
been understood as a motion characterized by the correlation between particle positions 
and momenta, of a dynamic origin. Early attempts to identify signatures of collective phe- 
nomena in single-particle observables were, generally, plagued by ambiguities. Thus it was 
anticipated [UJ|J that nuclear matter gets stopped in the very central collisions, compressed 
and heated, and subsequently expands primarily in the directions perpendicular to the beam 
axis. The expansion was suggested @ to be the cause of observed differences in the trans- 
verse slopes of the inclusive pion and proton spectra. However, it was further realized that 
the production of pions in the reactions mainly proceeds through the A resonance excitation 
and decay, and thus the pion slopes could reflect the A-decay kinematics ||. A sideward 
peaking of the proton angular distributions in an asymmetric reaction gave a hint of the 
expected collective behavior ||. As eventually, though, the neutron distributions had shown 
no such peaking 0, one was led to the conclusion that the peaking of the proton distribu- 
tions was a Coulomb effect ||. Stock et al. || noted that the measured pion multiplicities 
in central collisions, lower than predicted within the cascade or fireball models, could be 
explained by assuming that a part of the available energy was used up for compressing nu- 
clear matter. Subsequent dynamic calculations of collisions £[0| have demonstrated, though, 
a weak sensitivity of the pion yields to the compressibility of matter. 



Analyses of the 47r data brought fairly unambiguous evidence |TlJ of what appeared to be 
a different type of collective motion than first envisioned in the head-on collisions. This was 
a sideward deflection in the reaction plane of the fragments moving forward and backward in 
the center of mass in semicentral collisions. The deflection was quantified with a flow angle 
or with an average momentum in the reaction plane at a given rapidity. At midrapidity, 
at high beam energies, a preference was found for fragments to be emitted out of the reaction 
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plane |12,|l3|. It has not been obvious whether this preference was of a dynamic origin or 



due to shadowing. At low energies the emission in the reaction plane was found to be 



enhanced [TM 



Recent years brought new important results on the motion in central reactions. Specifi- 
cally, Schussler et al. [ T^IB \ have investigated the mean cm. energies of fragments emitted 



from the central reactions with Ag in emulsion. For the incident 36 Ar ions at 65 ± 15 MeV, 
and the incident 94 Kr ions at 95 ± 25 MeV, they observed a rise of these energies with frag- 
ment charge, faster than anticipated on the basis of Coulomb considerations alone. For the 
two reactions the rise was consistent with a collective expansion characterized by the energy 
of ~3 MeV/nucleon and > 7 MeV/nucleon, respectively. No rapid rise was observed for 
the incident 16 ions at 210 ± 10 MeV. The FOPI collaboration at GSI has further anal- 
ysed the spectra and mean energies of fragments with different charge emitted from central 
Au + Au collisions at beam energies in the range (150-400) MeV/nucleon |I7| , [l~8|j . They 
have found that the spectra of fragments with large Z were too broad and, respectively, the 
mean energies too high to be accounted for by a model of purely statistical decay. Within 
preliminary analysis [|17|| , data description improved once an assumption of the collective 



expansion at the time of decay was introduced, of an energy ~18 MeV/nucleon at the beam 
energy 150 MeV/nucleon, and of an energy ~51 MeV/nucleon at 400 MeV/nucleon, re- 
spectively. In a later more cautious investigation the collective energy was bounded 
from below with ~ 10 MeV/nucleon in the 150 MeV/nucleon reaction. The fragment yields 
were separately analysed to extract entropy [19[ within the so-called Quantum Statistical 
Model (QSM). Low entropy values were extracted consistent with low temperatures at rea- 
sonable values of freeze-out density. (These values were lower than those extracted in the 
past from light-fragment yields POl). The collective expansion could alleviate but not elim- 
inate discrepancies between the temperatures and energy available in the center of mass. 
Further analysis of spectra for collective motion in a central 100 MeV/nucleon Au + Au 
reaction was done in pi]. 



On the theoretical side, emission from central collisions has been analysed within a trans- 



port model with a dynamic light-fragment production (A < 3) f22| . Within this model the 
compression and excitation of nuclear matter in central collisions was observed, followed 
by an expansion. The mean energies of emitted fragments were observed to rise with the 
fragment mass, and a general consistency has been established between the rise with the 
mass and the collective energy determined from the dynamics. The calculated energies par- 
tially agreed and partially disagreed with the Plastic Ball data f23|. A degree of agreement 



of the calculations was found ||24fl , on the other hand, with the preliminary FOPI data on 
light fragments. Within the calculations ||22|| , the expansion affected, besides energies, the 
fragment angular-distributions. A midrapidity out-of-the-reaction-plane peaking at interme- 



diate impact parameters in a heavy system, in a quantitative agreement with the data (13 
(see also Refs. [^,^6[]), was identified as a continuation, with the rise of impact parameter, 
of a 90° peaking in the polar angle in head-on reactions. 



Spurred by the recent experimental evidence [15-|18|,|2T| and some forthcoming re- 



sults p7|,|2^j28[] , this study extends that carried out in Gil. Very preliminary results have 



been presented in [29||. Conditions under which a compression occurs in energetic reactions 
of heavy nuclei, and a collective expansion develops, are analyzed in detail, together with 
their effects on emitted light baryons and pions. Relations are established which approxi- 
mately govern the behavior of density and entropy in the compressed region in reactions, 
with beam energy and impact parameter. Anisotropies in particle emissions in the polar an- 
gle in head-on reactions and in the azimuthal angle in semicentral reactions, are explained in 
terms of a different starting time for expansions in different directions. Besides the effects of 
collective expansion on baryon and pion spectra, effects on particle yields are demonstrated. 
The effects on pions turn out to be dramatic, contrary to naive expectations. The role of 
expansion in the formation of hollow structures at lower energies p0|- p2| , |34| is elucidated. 



The analysis renders an understanding of both recent and older data. 

The outline of the paper is as follows. Section [TT] describes the model used for reac- 
tion analysis. Section |TJ is devoted to the dynamics and some observables from symmetric 
head-on reactions of heavy-nuclei. The gained experience is next utilized in the study of 



semi-central reactions in Sec. |IV|. Section |V| is devoted to the effects of collective expan- 



sion on pion yields, effects of the expansion on spectra of relativistic particles in general, 
and to shadowing of pions in semi-central reactions. Results on entropy and on collective 
expansion energy from measurements and simulations are discussed in Sec. [VJ. Section |V11| 
discusses the sensitivity of collective energy and its components to nuclear compressibility 
and elementary cross sections, and further the mechanism of formation of exotic structures 



at low energies. The results are summarized Sec. [VIII . 



II. TRANSPORT MODEL 

Reaction simulations are carried out within a transport model with explicit nucleon, 
deuteron, A = 3 cluster, pion, and delta and N* degrees of freedom |35| , |22[| . The phase- 



space occupation functions fx of stable particles satisfy transport equations following from 
the nonequilibrium many-body theory in the quasiparticle limit, 

at opx or or op x 

In the above E x is the single-particle energy, and K x and JC X are, respectively, the pro- 
duction and absorption rates of the particle X. The upper and lower signs in the statistical 
factor on the r.h.s are for fermions and bosons, respectively. 



The deuterons within the model are formed in the interaction of three nucleons 35,37.38 



The formation rate @,[3lJ contributing to the deuteron production rate has a form analogous 



to that of the rate associated with elastic scattering, 

rrid ^ f dp m x dp\ m x dp' 2 m x dp' 
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The dots indicate other terms in the production rate. The factor \Ai\ 2 stands for the 
matrix element squared summed over the final and averaged over the initial spin directions. 



A term in the deuteron absorption rate accounting for the deuteron breakup into nucleons 
has an analogous form to the term for the formation (0), with statistical factors for the 
initial and final states interchanged. Due to the microscopic reversibility the processes of 
deuteron formation and breakup share the matrix elements squared. The averaged elements 
squared are then related with the ratio of spin-degeneracies, 

3: 



l-M-Npn^Ndl 2 — -^{■M-Nd^Npnl 2 - (3) 

In principle, data on deuteron breakup could be used in describing the formation as the 
differential cross section for breakup is proportional to the matrix element squared, 
1 m N 



da Nd ^ Npn = —^r-f — r\M Nd ^ Npn \ 2 27rS(m d + £jv(pi) - E N (p[) - E p (p' 2 ) - E n (p' 3 )) 

In practise, due to the dimensionality of three-body space, the differential cross section is 
known only within a limited range. In the calculations, a modified impulse approximation 
is thus used 



\M pd ^ ppn \ 2 ~ f{ Kp'MI 2 \M„ n \ 2 + Kp'MI 2 \m 



pn-^pn | 



+ l(P3l^>| 2 |M^ PP | 2 }, (5) 
where <fi d is the deuteron wave-function in momentum space, normalized so that 

The three terms in (|5]) correspond to the three different nucleons being spectators. 
The NN matrix elements in are proportional to the NN cross sections. Finally, the overall 
normalization factor F in ([5]) is adjusted so that the measured total deuteron breakup cross 
section is reproduced as a function of bombarding energy [|S5|| . 

Tritons and helions in the transport model are produced in the interactions involving 
four nucleons p6| , |22| . The production of a composite is generally suppressed if the average 



nucleon occupation over a volume in momentum space, corresponding to the composite wave- 
function, exceeds a phenomenological cutoff-value of 0.30. At low densities, in the limit of 



many interactions, the set of equations for the composites and nucleons yields the required 



law of mass action. A coalescence scaling [39] for the composite spectra is expected to 



emerge from the model based on the equations, when the processes of composite production 
do not vary substantially with space and time. 

The pions within the model are produced in two steps [fHJ. First a A or iV* resonance 
is produced in an NN collision, and then the resonance decays into a nucleon and pion. 
The cross section for the resonance production appears to nearly exhaust the inelastic NN 
cross section in the energy range of interest ||41|| . The resonance occupation functions satisfy 



transport equations that follow from the nonequilibrium many-body theory in the adiabatic 



limit 35 



Wt {fx Ax) + 9pf 'di {fx Ax) ~ "dT'dpT ifx Ax ) = 1C x^T fx) A x - K> x f x Ax, (6) 
where Ax is the spectral function describing the spread of the resonance in mass, 



T x 

(m-m x ) 2 + \T x 



and Tx is the resonance width which generally depends on m. The derivatives on the 
l.h.s. of Eq. (Q) are taken at a constant m — mx- The absorption of pions is described 
with a sequence of inverse processes to those in the production. The resonance formation 
dominates the ttN cross section at the energies of interest, with the cross section for resonance 



X of the form 42 



71 9x r , , Q , 

p z 2 

with p - cm. momentum and gx - resonance spin-degeneracy. Microscopic reversibility 
and the assumption of a weak dependence of the interaction matrix-element squared on the 
resonance mass, yield a relation between the cross sections for resonance absorption and 
production p5| , ^3| , 

12 mp 2 N , N „ <j n < n »^xn /ri N 

axN^N'N" - — -= E _ mN j^rz — , I9j 

1 + Ojv'jv" gx Pxn J amAximjpxN 



where E is total energy. A variant p4j of the relation relies on the assumption of a rapid 



variation of the matrix element, for m close to E — m x . In the calculations preceeding ||35|| , 
a detailed balance relation between the resonance production and absorption cross-sections 
used was such as for stable particles, leading to an excessive number of pions in equilibrium. 
In the pion-nucleus interactions the true pion absorption cross section was underestimated. 
Further violating reversibility in the calculations, in (||) the momentum prefactor was re- 
placed by a constant, affecting the shape of pion spectra. 

The single-particle energies in Eqs. (|l]) and @ are parametrized in terms of an optical 
potential, with a local part handled as a scalar potential, and a weaker nonlocal part added 
on to the particle energies in the system cm., 

E x = {p 2 + m 2 ) 2 + q x $ + A x U? , m = m + A x U c + t 3X U T , (10) 

where q x is particle charge, $ - Coulomb potential, A x - particle mass number, m - mass 
in free space, and t 3X - isospin component. The dependence of the potentials U on scalar 
densities is chosen such as in nonrelativistic calculations, 

uC = - a 7 + b (fi- u ° = ~ dv2 (?) ' " T = cP > 

where 

Ps = ^ZA x p sX , p J = J2t 3 xpsX, (12) 

x x 

with 

f dp m xo f dp fdmm xo 

for the stable particles and resonances, respectively. The neutron and proton density profiles 
for initial nuclei are determined by solving a differential equation for baryon density that 
follows from the Thomas-Fermi (TF) equations in the nuclear frames, 

E p (r,p Fp (r)) = n P , E n (r,p Fn (r)) = \i n . (14) 
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The chemical potentials /i p and /i„ are adjusted till the values of total proton and neutron 
numbers for nuclei are reproduced. It appears convenient to solve differential equations set 
up for the proton and neutron numbers, Coulomb potential, and net energy in parallel to 
the equation for baryon density 

The normal density is taken equal to po — 0.160 fm~ 3 . The description of charge densities 
from electron scattering ||5| is better with this value than with po — 0.145 fm -3 ||22j| . As to 
the other parameters in (0), a = 357 MeV, b = 304 MeV, and v = | correspond to the 
soft equation of state (K ~ 200 MeV), and a = 123 MeV, b = 70.1 MeV, and v = 2 to the 
stiff equation (K ~ 375 MeV); c = 92 MeV, d = 22 A MeVfm 2 |l,|47j. Sample densities 
for the second parameter set are shown in Fig. [I]. Unless otherwise indicated other results 
in the paper are obtained using the second set. The effects of the momentum dependent 
potentials on the expansion have not been explored as of yet, due to the involved numerical 
difficulties with energy conservation and a specific interest in the differences between mean 
energies and spectra of different emitted fragments. 

Occupation functions in the phase space for stable particles, and in the phase space and 
mass for resonances, are represented in the calculation with test particles. The drift terms 
on the l.h.s. of Eqs. ([T]) and (g) are integrated by requiring that the test-particle positions 
and momenta satisfy the Hamilton's equations. A lattice hamiltonian method of Lenk and 
Pandharipande f46[ is used with the scalar density at lattice sites evaluated with a form 
factor 

S x (r,p) = A * m *° g(x)g(y)g(z), (15) 

l\test{2t)° hj X 

where g(q) = 1 for \q\ < 0.5£, g(q) = 1.5 — \q\/£ for 0.5£ < \q\ < 1.5£, g(q) = for 
|g| > 1.5£. Further, N test is the number of test particles per particle, and £ is lattice 
spacing. The particular form of g allows for a good total energy conservation (to within 
0.25 MeV/nucleon in a 100 MeV/nucleon collision), and a good momentum conservation for 
asymmetric systems, with no considerable computational effort. It should be mentioned that 
an initialization of the system according to the densities from the TF equations is essential for 
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an "elastic" algorithm for the drift terms. For other, however sensible, initializations, isolated 
nuclei develop strong breathing oscillations that can deplete central density to a quarter of 
normal value. The Coulomb potential is found by solving the Poisson equation using a novel 
relaxation method . Integrals in the collision rates in Eqs. ([[]) and (|6]) are computed using 
the Monte-Carlo method (see also Ref. Entropy is calculated by integrating |J5 

the rate of entropy variation due to collisions. 



III. HEAD-ON COLLISIONS OF HEAVY NUCLEI 

Emphasis is going to be on symmetric reactions of heavy nuclei. For symmetric systems 
the experimental identification of specific collective effects may be less ambiguous than for 
asymmetric systems, due to the well defined center of mass and the lack of shadowing at 
lowest impact parameters. 



A. Shock Fronts 

The dynamics of the central high-energy reactions can be broken down into several well 
identified stages. A 400 MeV/nucleon Au + Au system at b = will be used to illustrate 
specific points. Contour plots of baryon density at different times for that system are 
shown in the left column of Fig. |2|. Initially in the energetic reactions the nuclear densities 
just interpenetrate, with target and projectile nucleons separated in the momentum space. 
As the time progresses, in the overlap region NN collisions begin to thermalize the matter, 
making the momentum distribution centered at zero momentum in the c.m.s. Thereafter, 
the features of the matter at the center of a system stabilize to a degree. As further the 
region with the excited matter grows in size, in a heavy system interface regions can be 
identified in-between the excited and normal matter where such parameters as baryon and 
entropy density change rapidly, see Figs. [2] and []. In the head-on collisions the interfaces 
are perpendicular to the beam axis. With the normal matter diving into the region with 
excited matter at a speed well in excess of the speed of sound in normal matter (for the first 
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sound the speed is c s = yjK/9rri]y, i.e. 0.15c and 0.21c in the case of soft and stiff equation 
of state, respectively), at the high beam energies, the interfaces are recognized as shock 
fronts. In the hydrodynamic limit a discontinuity in the velocity in the initial state, such as 
in-between the projectile and target nucleons at b = 0, generally breaks into two shock fronts 
travelling in opposite directions [50]. Figure § displays the nuclear-matter parameters within 



the hydrodynamic limit, behind a shock front bringing the matter to a halt, as a function 



of beam energy, obtained by solving the Rankine-Hugoniot (RH) equation [50 

(ei + PO 2 e% (e 1 + P 1 , e 



, ^R^ + i =o. (16) 
Pi Po V Pi PoJ 



The parameters with subscript 1 in ( jig) refer to the shocked matter and e and P are the 
energy density (including mass) and pressure, respectively. Figure |] displays further the 
maximum densities from the 6 = Au + Au simulations and these densities appear to 
follow the general trend predicted by (TO). The excited matter in-between shock fronts in 
the simulations is exposed to the vacuum in the sideward directions and, respectively, begins 
to expand into these directions. The features of this process at high energies and at b = 
may be understood at a qualitative level in terms of the self-similar cylindrically-symmetric 
hydrodynamic expansion of an ideal gas. 

B. Expansion 



In the self-similar expansion [51,52], the velocity is proportional to a distance from the 
symmetry axis, v = J-{t) r, which roughly appears to hold in the direction of 90° in the 
Au + Au collision at 400 MeV/nucleon, cf. Fig. Under the simplifying assumptions of 
an isentropic flow and an ideal-gas equation of state, yielding Pj p 1 = const with 7 - the 
heat ratio, the solution to nonrelativistic hydrodynamic equations for the density becomes 

r 2 (o) ( r 2 \ 1/(7 - x) 

P^)= P (0,0)^[l- m ) . (17) 
The radius R of the density distribution satisfies 
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~ _ 2^(0,0) (R(0)\*~' 

m - {Wi) • (18) 

with (? s = jT/m - the speed of sound squared. The coefficient JF is related to R with 
J 7 if) = R(t)/R(t). Equation QT8D may be partially integrated to yield, for R{0) = 0, 

27-l N 



R(t) = c s (0, 0) 



7(7-i)l Uw; 



(19) 



Following Eqs. (|l8j) and (|i~9l) , the coefficient T is expected to rise, first linearly, for the time 
t ~ R/c s needed for a signal to propagate inward through the radius, and it is expected to 
fall off later as the inverse of time, reflecting the quadratic initially and linear asymptoti- 
cally dependence of R on time. Indeed, in the high-energy simulations, the velocity slope 
dv/dr in the transverse direction, cf. Fig. |3], behaves approximately as described. Within 
the nonrelativistic ideal-gas approximation, the speed of sound in the central region of the 
collision at 400 MeV/nucleon is assessed, with 7 = 5/3, at c s ~ 0.4c. The radius of the 
gold nucleus is ~ 7 fm, and thus the slope might be expected to rise for a time ~ 18 fm/c. 
The maximum in the simulations is reached at t ~ 27 fm/c which appears consistent with 
the simple estimate as initially some time is spent equilibrating and compressing matter, 
cf. Fig. [5]. The value at maximum dv/dr m 0.040 c/fm is not far from the one obtained 
by integrating (|i~9l), T ~ 0.6 c s /R ~ 0.036 c/fm. At the given instant, according to (|l~9|). 



the transverse matter radius should be about 50% in excess of the original which approxi- 
mately is the case in the simulation, cf. Figs. [| and [3]. After the time R/c s , the density in 
the central region of the reaction should decrease as t~ 2 , following Eqs. (|17D and (|T9|). 

While the expansion develops in transverse directions, the normal matter continues to 
enter the excited region along the beam axis for a time ~ (2R/vo)(pi — po)/pi ~ 16 fm/c in 
the Au + Au collision. Here Vq ~ 0.42 c is initial velocity in the center of mass and p\ is the 
density in the shocked region. After the shock fronts reach the vacuum, an expansion along 
the beam axis sets in that requires a time ~ (2i?/c s )(po/pi) ~ 18 fm/c for a full development. 
By then, however, the matter is already quite decompressed with particles having moved 
out in the transverse directions. At the later stages, with the expansion becoming three- 
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dimensional, the density in the central region can be expected to drop as t~ 3 . 

While the discussed and few other dynamic features of b = high-energy reactions 
appear to be consistent with a hydrodynamic behavior of the matter, there are also definite 
differences which cannot be classified as effects of the mean-field potential or relativity. Thus, 
the baryon density at the center of the system, Figs. |3|, and |5], falls short of the value 
of 2.5 po expected at 400 MeV/nucleon on the basis of the RH equation. The entropy per 
baryon never quite reaches values expected for an equilibrated system at the center. Typical 
combinations of density and mean nucleon kinetic energy such as p ~ 2.1 po and Eki n ~ 
91 MeV, should result in Sb/A ~ 2.8, compare Fig. [| For reference, the nonrelativistic 
formula for the entropy, with statistics ignored, is 



Sb ~ - . I P ( 37rfr 



^ = 2.5 - log | ^ T I I. (20) 



■2 \3/2N 



The shock fronts do not fully develop |)3] and/or the system never gets completely ther- 
malized at the center (see also ||54|| ). In a system expanding hydrodynamically into vacuum, 
local temperatures would have dropped to zero and all kinetic energy would have turned 
into the collective energy of expansion J5l],[52|. For that, collisions would need to continue 



down to appropriately low densities. The mean density in the vicinity of last collisions of 
particles emitted from head-on Au + Au reactions, is shown as a function of the beam energy 
in Fig. H At 400 MeV/nucleon the freeze-out density is equal to about 0.2p - The bottom 
panel in Fig. |5] shows the evolution of the participant collective transverse energy per nucleon 
in the 400 MeV/nucleon reaction, which is calculated from the perpendicular components 
of a local hydrodynamic velocity v in the vicinity of last collisions, 

Eiu/A = £ / drdpfxmx^) 2 / drdpfc . (21) 
x J I x J 

The collective transverse energy per nucleon freezes out in Fig. |5| at a value of 42 MeV. 
For nucleons this constitutes about 50% of the kinetic energy in transverse degrees of freedom 
at the freeze-out. (A bit larger fraction is found if relativistic effects are taken into account.) 
The nuclear systems are, indeed, too small for the entire energy to turn into the collective 
energy. Simple estimates of mean freeze-out density based on 
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where r = (pctnnv and t>* is relative thermal velocity, that utilize a nuclear density 
p{t) oc t - " where v — (2 — 3), and rely on the assumptions of either adiabatic or isothermal 
expansion yield (p) / pi oc (<jnnRpi)~ k , where p\ is the density at the start of expansion and 
k ranges from 2/3 to 2. Following these estimates, the energy dependence of the freeze- 
out density seen in Fig. ^, above 100 MeV/nucleon, may be attributed to the variation 
of elementary cross sections of particles interacting locally with reduced relative velocities. 
If cross sections are made energy-independent in simulations, indeed little change in (p) is 
observed towards high energies. 

One of the effects of transverse expansion on observables from very central collisions 
in the simulations, is the peaking of particle polar-angle distributions at 90° in the c.m.s. 
This peaking is more pronounced in the distributions of clusters which are more affected 
by the hydrodynamic motion than are nucleons | 55| , see e.g. Fig. |7]. With the variation 



of beam energy, the peaking is observed to maximize at about 150 MeV/nucleon. Be- 
low 40 MeV/nucleon the emission turns isotropic. The shock width at low energies becomes 



comparable to the nuclear size |53[ and thus a system equilibrates as a whole. For polar- 



angle distributions from the simulations of 250 MeV/nucleon reactions see Ref. p2[ . At very 
high energies (and/or for light systems like Ca + Ca) the effects of expansion in the simula- 
tions compete with the transparency effects in the corona. The proton distributions become 
forward-backward peaked in the c.m.s., at beam energies above 1.5 GeV/nucleon, and the 
light-cluster distributions above 1.9 GeV/nucleon in the Au + Au reactions. As seen in 
Fig. |7|, the pion polar-angle distribution peaks at 9 = already at 1 GeV/nucleon. For the 
purpose of demonstrating that, indeed, the transparency plays a role, the figure shows, in ad- 
dition to the distribution of all pions, the distribution of only those pions which were emitted 
within the first 14 fm/c. The latter pions constitute about 1/3 of all, and the anisotropy 
of their polar-angle distribution is nearly 3 times larger than the anisotropy for all. (Possi- 
ble effects of geometry in a surface emission in shaping the pion angular distribution may 
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be eliminated, as the angular distributions at lower energies are peaked at 90°.) Despite 
the transparency effects in the corona, the maximum densities at the center of the Au + 
Au system are not farther away from the RH expectations at the highest of the studied 
energies than at lower,[] cf. Fig. The peaking of particle distributions at 90° at b = 0, 
including that for pions at lower beam energies, eventually disappears with an increase of 
impact parameter f22[. Experimentally, the most central collisions can be selected by using 



the particle transverse energies or symmetry around the beam axis ||18|| . In order to avoid 
a bias in the angular distributions, the analysed particles should be excluded from the trigger 
condition (compare f56fl). By forcing the condition to a limit the extraction of information 
from most central collisions could be reduced to the issue of collecting an adequate event 
statistics. 

Figure H shows the mean transverse and longitudinal kinetic-energy components p2]j , 



E kin = E-m = Jm 2 +p 2 - m = y - ' + y - ' = E ± + E\ (22) 
km v y m+E m+E v; 

of particles emitted from the 6 = Au + Au collisions, as a function of the bombarding 
energy. The scale for the transverse component is shifted downward by log 2 relative to 
the scale for the longitudinal component, in order to allow for a direct visual comparison 
of the components per degree of freedom at the same bombarding energy. It is seen that, 
per degree of freedom, the transverse baryon energies are higher than longitudinal in a vast 
range of bombarding energies. This is a reflection of higher, per degree of freedom, collective 



transverse energy than longitudinal, at freeze-out f2~2 |. The components of the collective 



energy per nucleon, calculated as in Eq. (^1|), are shown in Fig. |9|. If particle momentum 
distributions are approximated with anisotropic gaussians and nonrelativistic kinematics is 
used, then the particle polar-angle distributions become 



1 In a comparison in \22\ over a more narrow range of energies, the pions and baryon resonances 



were omitted in the RH prediction and some nonrelativistic approximation was made. Besides, the 
predictions there were for another equation of state. 
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47T dN (1 + X) 1/2 m , 
N dtt (1 + xcos 2 ^) 3 / 2 ' 1 ' 

where x = ((-^ ± )/2 — (E"))/(E^). Result (p3j) for tritons emitted from a 1 GeV/nucleon 
Au + Au collision is shown by a solid line in Fig. [7[ In general, it is found that 



Eq. (23) describes the distributions of emitted particles from high-energy collisions in a semi- 



quantitative way, even the pion distributions. The 90°-to-0° ratios, 



'dN\ /( dN \ 



are quite well described within the approximation, Rq « (1 + x) 3 ^ 2 = ((E^) /2(E^)) 3 ^ 2 . 
Within the gaussian approximation, the mean particle energy at 90° is equal to 1.5 x 
(E- 1 ). The mean energies at 90° in simulations are slightly lower. 

IV. SEMICENTRAL COLLISIONS 

The attention is now turned to the dynamics of semicentral collisions, with emphasis put 
on the similarities and differences compared to the dynamics at b = 0. 

A. Hydrodynamic Discontinuities 

The initial discontinuity between the projectile and target velocities breaks at a finite b 
in the ideal-fluid limit into two shock fronts propagating into the projectile and target, 
and a, weak tangential discontinuity in-between |EHBD[. Unlike across a shock front, the nor- 



mal velocity component, pressure and, in the case at hand, entropy and density are con- 
tinuous across the tangential discontinuity, and the tangential velocity component is not. 
The discontinuities start out in collisions all at the same angle a relative to the beam axis, 
cf. Figs. |ig and g 



b 1 
cos ct ~ — (24) 

2R Jl + (E beam /2m N )(l + byAR^ 
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where the beam energy is in MeV/nucleon, and the Lorentz contraction in the c.m.s. is 
taken into account (v = (1 + 2m N / E beam )~ l l 2 , 7 = J 1 + E beam /2m N ). Developing in- 



clined shocks can be seen in the t = 10 fm/c panels of Fig. |2| for a finite b. Signs of 
a tangential discontinuity for a finite b are visible in the densities at later times in the figure. 



Figure [TT| shows finally the variation of different quantities along a normal to discontinuities 
in a collision at b = 7 fm. Note that the center of the weak discontinuity is at r = where v t 
vanishes, while the center of a shock is at r ~ 3.1 fm where v n reaches 1/2 of the value for 
the intact matter. 

In the ideal-fluid limit, the tangential component of a velocity would be unaltered by 
a shock, and thus the matter behind the shock in collisions would be expected to move at 
a speed v\ = v cos a and angle a relative to the beam axis. The change in the thermo- 
dynamic parameters should be then such as for a shock perpendicular to the beam axis, 
at b = and a beam energy (called further an equivalent energy) 

Keam = E beam S1H 2 a, (25) 

where relativity has been thoroughly accounted for. The ideal-fluid dynamics predicts then 
that the maximum density reached early on in collisions in the cm. at a finite b, is related 
to the density at b = at a lower beam energy, by 



pi(6, E beam ) = 7ipi(0, E' bean 



V 1/2 „ fn „ (l+e)(l-/3y 



Pi 0, E beam - 



V (i + e) (i + e(i + p))j r± v eum i + 6(i+/3) ; 

pi(0, E beam (l - 0)) = Pl (0, E beam (l - 6 2 /4i? 2 )), (26) 



where e = E beam /2mN and (3 = b 2 /4R 2 . The approximation following the equalities in ( |26|) 
is nonrelativistic (e 1). The cm. maximum density and corresponding rest-frame pa- 
rameters like entropy behind the shock fronts are expected to change slowly according 
to (|25"D and ([24]) , only quadratically, with the rise of b. Indeed, a weak variation behind 



the shock fronts is found in the simulations. The validity of the relativistic relation ( 26 ) 



is tested in the top panel of Fig. |4] that displays, in particular, the maximum density from 
1 GeV/nucleon Au + Au collisions at different b, divided by 71, vs E' beam from (|26|). The re- 
lation appears to work up to b ~ 11 fm at this and other beam energies. Further predictions 
based on the ideal-fluid dynamics need not work as well. 

In the simulations all the discontinuities have widths that are quite comparable to the 
size of a system, see Fig. [0]. The tails of the weak and strong discontinuities overlap 
and in consequence the velocity in-between does not, in fact, reach values such as if the 
discontinuities were separated. For example, at low b, along the axis passing through the 
center, the magnitude of the transverse velocity component in the reaction plane reaches 
at most 40% of (v sin a cos a) expected in-between separated discontinuities. As a function 
of time, the width of any shock front (of the order of a mean-free path for a strong shock) 
remains constant |53|j50[| , provided that the conditions on the two sides of the front remain 



the same. The width of a tangential discontinuity, on the other hand, increases with time. 
That this must be the case can be seen directly by adopting the hydrodynamic Navier-Stokes 
equations [0] to a region of the discontinuity, where they yield a diffusion equation for the 
tangential component of velocity (or for vorticity u = V x v). The diffusion coefficient is 
equal to r]/{m^p), where 77 is the shear viscosity coefficient (proportional, in general, to the 
mean-free path). Solving this diffusion equation yields the tangential velocity component as 
a function of time t and distance r in a direction perpendicular to the discontinuity, in the 
form 



(27) 

where is an asymptotic value of the velocity component and erf is error function. Accord- 



ing to Eq. (p?D the half- width of a weak discontinuity should rise with time as 5 = 2yjr}t/m^p. 
For the situation in Fig. O], with 77 ~ 50 MeV/fm 2 c |[53|| , this yields 5 ~ 3.2 fm which agrees 
well with what can be deduced from the simulation in the figure. Generally, during the 
time span lasting from the initial equilibration in simulations till the development of the 
expansion at intermediate and low impact parameters, the tangential- velocity profile at the 
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center, is reasonably well described by Eq. fl27|), despite of the proximity of shocks. The root 
in ( p?D alleviates uncertainties associated with the choice of density and a reference time. 
While 5 increases as t 1//2 , the separation between shocks increases as t, and thus, if nu- 
clear systems were large enough and/or viscosity low, the shocks in collisions could separate 
well from the weak discontinuity, before reaching the vacuum. For large Reynolds numbers 



1Z = rriNpvR/ri the weak discontinuity at the center could break into a turbulence J50 



In the system studied TZ ~ 9, while typical critical values for turbulence are larger |50 
No sign of any turbulence is seen in the simulations. 

In the course of a reaction, the section of the weak discontinuity turns from one given by 
a straight line to one reminding a hyperbolic tangent, and on the average the discontinuity 
rotates, decreasing an angle relative to the beam axis. That is because the matter moving 
forward on one side of a reaction plane, away from the center of the system, emerges facing 
little or no target matter moving the other way. Symmetrically, the target matter on the 
other side emerges facing little or no projectile matter. The shock fronts subsonic relative to 
the matter behind rotate as well. The fronts are generally strong in the vicinity of an axis 
passing through the center of a system in the direction of the beam; there from the shock in 
projectile matter continues as strong towards the target side of the reaction plane, and the 
shock in target matter as strong towards the projectile side of the reaction plane. With the 
excited region within the geometric overlap of projectile and target fed by incoming matter, 
the density is maintained behind the shocks during the shock propagation. 

The shock in projectile matter is weakest in the vicinity of its edge on the projectile 
side of the reaction plane. There the front curves, making a low angle relative to the beam 
axis, the change in the velocity component normal to the front is low, and correspondingly 
the equivalent energy for the front is low, compare (^). Symmetrically, the front moving 
into target matter is weakest in the vicinity of its edge on the target side of the reaction 
plane. The distinction between the strong and weak portions of the fronts is generally 
pronounced at high beam energies but less at low. The weak shocks propagate with only 
slightly supersonic velocities relative to the intact matter. These shocks, in what becomes 
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spectator matter at high energies, might remind those away from supersonic bodies moving 
through a medium. The matter little disturbed by the weak shocks can move out as an entity 
largely in its original direction in the cm., see the t = 35 fm/c panel in Fig. ^| for 6 = 7 fm. 
(The spectator matter separates out also in the reaction at b = 3 fm at a later time than 
shown in the figure; the shapes for the matter are far from spherical.) 

If the hydrodynamic discontinuities in the region of geometric overlap between projectile 
and target, were more separated in the simulations, the matter would have moved more 
to the sides within the reaction plane behind the shocks, strengthening the shock edges. 
The spectator matter would not have been then as distinct; see the results in j58]|. Weakness 



of the sideward flow in collisions [p9| , as compared to the ideal-fluid predictions, and the 
emergence of spectator matter are thus very directly related. 



B. Expansion 

From the shock formation on, an expansion proceeds to develop in-between the fronts. 
At low b the expansion in the directions perpendicular to the reaction plane has features sim- 
ilar to the expansion in the directions perpendicular to the beam axis at b = 0. The variation 
of baryon density, entropy and radial velocity at different time instants, with the distance 
from the center of the system along a line defined by 9 = 90° and <fi = 90°, where <p is the az- 
imuthal angle relative to the reaction plane, is quite close to the variation given by the dashed 
lines in Fig. |3|. The top panel in Fig. [12] shows the dependence on the impact parameter of 
the mean particle kinetic energy components out of the reaction plane E y = (p y ) 2 / (m + E) , 
where p y is the momentum component out of the reaction plane, at cm. rapidity y = 0. It is 
seen that the mean energy components depend roughly quadratically on b for low values of b 
(i.e. in such a manner as the maximum density in the collisions). Center panel in the figure 
shows the energies E y averaged over whole momentum space. These averages are equal to 
half of (E- 1 ) at b = 0, and they decrease faster with an increase of b than the averages at 
y = 0, because of larger spectator contributions. Either type of average energies may be 
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determined experimentally, correcting for fluctuations of the reaction plane [^2| , p0 , p6 . 

With regard to the directions in the vicinity of the reaction plane, the initial velocities 
tangential to the weak discontinuity behind shock fronts (although low compared to those in 
the ideal-fluid limit) add on to the velocities developed in expansion, and affect the sideward- 
flow observables. The deflection of particles moving forward and backward in the cm. is 
commonly quantified in terms of the mean momentum in the reaction plane as a function 
of the rapidity [p6[l . Figure shows the mean momentum divided by the mass, (p x {y)) /m, 



for particles emitted from Au + Au reactions at 1 GeV/nucleon at two impact parameters. 
For the purpose of minimizing the effects of spectator particles at high energy, the slope of 
the mean momentum at midrapidity is further used in analyses The slope 

d(p x /m) 



dy 



(28) 



is shown in Fig. [TJ] for 400 M eV/nucleon Au + Au collisions, as a function of the impact 
parameter, both for the calculations and data. This slope generally increases with particle 
mass. If the emission pattern at low h and low cm. rapidities differed from that at b = only 
by a rotation of the plane of preferred emission, from the angle of 90° relative to the beam 
axis, to an effective angle a, then the slope F at midrapidity should be approximately [ |T2"| , p2] 

,_, y sin a cos a . . 

F « f- — , 29 

with x obtained from the energy components at b = 0, cf. (|23|) . In obtaining (|29|) non- 
relativistic kinematics was used and momentum distribution in the vicinity of midrapidity 
was approximated in a gaussian form. The maximum value of F from (p9|) as a function of 
a equals to xl (2\A + x)- F° r protons from 400 MeV/nucleon collisions this is 0.22 (using 
energy values from Fig. §). However, the proton F in Fig. [14] reaches ~0.47 at intermedi- 
ate impact parameters, indicating the importance of the weak discontinuity. Note that the 
spectator contributions would tend to lower F. 

On taking spectator contributions into consideration, and on analyzing the behavior of 
the r.h.s. of Eq. ( P5Q with one finds that F provides a lower limit to the angle of 
inclination of the discontinuities within participant matter relative to the beam axis, 
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a > tan 1 F. 

The calculations in Fig. [14], utilizing parameters for the stiff equation of state and free 
NN cross-sections, actually somewhat overestimate data. Thorough discussions of the de- 
pendence of observables characterizing the sideward deflection on the microscopic input 



parameters in simulations can be found in Refs. |63| . It follows from Fig. |T^ that the highest 
multiplicity interval within the experiment corresponds to b< 2.5 fm. 



Either data or calculations may be analyzed in terms of a kinetic-energy tensor ||64j| . This 
tensor can be written in a relativistic form as 

S ij = y ij = y PlPi / 30 ) 

where the sum is over particles from a collision. Experimentally, the mean elements of 
the particle tensors (s y '), associated with the reaction plane, may be determined directly as 
outlined in [0 , irrespective of the fluctuations of an estimated reaction-plane about the true 



plane.0 The major axis of the mean particle-tensor points in a direction within the reaction 
plane, at an angle Of relative to the beam axis. This direction is intermediate between the 
direction of discontinuities within the participant matter and the direction of motion of the 
spectator matter. Eigenvalues associated with the axis pointing out of the reaction plane are 
identical with the mean energies (E y ). The minor axis within the reaction plane, denoted 
with symbol x', is generally at a considerable angle > 50° relative to the strong and weak 
discontinuities within the participant matter. As b — > this axis becomes identical with the 



2 Following considerations such as in Ref. optimal weights in a vector [56] for the estimation of 
the reaction-plane direction, Q = J2u w vPvi are w v ~ (p x ) I ({p ± ) 2 }, where averages refer to specific 
rapidities and particle types. A more general equation for the estimation of the azimuthal angle $ of 
the reaction plane follows from the requirement = ^ J2ui(Pu cos(4> u — &) — (p x )) 2 / {(p x — (p x )) 2 ) + 
sin — <£)) 2 / ((p y ) 2 })- The averages within the equations may be determined self-consistently 
in data analysis. 
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beam axis. Eigenvalue energies (E x ) are generally less than (E y ), see Fig. |T^. The ratios 
(E y )/{E X ) increase with particle mass. These ratios maximize at b = as a function of the 
impact parameter, where they coincide with {E- L )/2{E^). Compared e.g. to the variation of 
the flow angle Of (see also Ref. ||66||), the variation of the discussed energy ratios with b is, 



though, slow, cf. Fig. O and further Fig. 15 



Kinetic energy tensor, in its nonrelativistic form, was determined in the Au + Au reac- 



tions on an event-by-event basis in Ref. ||67|| . The mean ratio R\ was calculated of matrix 
elements associated with the direction out of the estimated reaction plane and the direction 
within the plane, perpendicular to main axis, as a function of participant proton multiplicity. 
Fluctuations of the estimated reaction plane about the true reaction plane generally reduce 
the possible differences between matrix elements in question ||56|| . (If not for fluctuations, 
R\ would have been identical with the energy ratio (E y )/ {E x ').) For the Au + Au reaction 
at 400 MeV/nucleon, the average azimuthal angle between the two planes has been as- 



sessed within the experiment, and given in Ref. |67j as function of proton multiplicity. This 
makes possible a practical comparison of present calculations with these data. Ratio of the 
eigenvalues from the tensor with matrix elements modified on account of the fluctuations!] 



is compared to the data |)7| in Fig. [15]. As in the case of Fig. [13] a general agreement is found 



although, quantitatively the calculations overestimate somewhat the data. The lowering of 
the ratio R\ at high multiplicities stems from an increase in fluctuations. With regard to the 
flow angle Of, the position of a maximum in the experimental event distribution in the an- 
gle p8| , for the highest multiplicity bin, appears to be in a rough accord with the conclusion 
on the impact parameter range, reached by examining F. 

Physical reason behind larger energy values out of the reaction plane in simulations, 



than in a direction within the reaction plane, such as in the case of Fig. [15], is the same 



3 The difference ((s xx ) — {s yy )) is reduced by a factor (cos 2(5$) and the element (s xz ) is reduced by 



(cos (5$), where <5<E> is azimuthal angle between the estimated and the true reaction plane [56 
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as behind the different energies at b = 0. The expansion of largely equilibrated matter 
in the direction of motion of shock fronts is delayed, compared to the direction out of the 
reaction plane, by the time necessary for the shock fronts to traverse the normal matter. 
As matter is already being decompressed by the expansion developing in directions parallel 
to the shock-front surface, the expansion in the direction of shock motion does not acquire 
same strength. The fact that data |)7|] (also fl2] , |60| ) show different energy values associated 
with different directions demonstrates that, indeed, the equilibration occurs within a part 
of the system separated out from the cold matter. With motion being generally supersonic 
at a high energy, the interfaces between the hot and cold matter must have the nature of 
shock fronts. The sideward flow stronger than anticipated on the basis of a rotation of the 
b = pattern, indicates the presence of a weak discontinuity (for the 6 = pattern an 
extrapolated energy ratio from finite-6 data such as in Fig. [15] could be used). 

For other beam energies than 400 MeV/nucleon, azimuthal deviations between the esti- 
mated and true reaction planes were not given in |J7|], and, correspondingly, the measured 
and calculated energy ratios of interest in the Au + Au reactions cannot be directly com- 
pared. Comparisons to the theory would have been eased if energies associated with the 
directions relative to the true reaction plane were established, as was discussed. Unless one 
deals with some artifact of a strong variation of the fluctuations with beam energy, a con- 
siderable discrepancy between data and calculations, with regard to the energy ratios, may 



be arising at low beam energies. Specifically, the experimental squeeze-out ratio R\ [67 



maximizes at 250 MeV/nucleon in central reactions and rapidly decreases from there on 
with a lowering of the beam energy. On the other hand the ratio of calculated energies 
(Ey)/{E X ) (which may be obtained directly from Fig. ||), maximizes at 150 MeV/nucleon 
and continues as significantly larger than 1 down to 50 MeV/nucleon. 
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V. PION PRODUCTION 



It is generally believed that, in comparison to baryons, pions are rather insensitive to 
the collective motion, due to their low mass and corresponding high thermal speeds. First, 
briefly, the possible effects of a collective expansion on energies and spectra of particles with 
different mass are investigated by considering a situation of the instantaneous freeze-out in 
the c.m.[] of a locally equilibrated system at a constant temperature. 

A. Collective Expansion and Mass 

With inclusion of particle statistics, phase-space distributions in a locally equilibrated 
system are of the form 

fx = ! 

expftvTpv -y. x )/T\ =F 1 
where the upper signs refer to Bose and the lower to Fermi statistics, respectively, u° = 
(7,7V), v is velocity, and fix is chemical potential. The limit of a classical statistics corre- 
sponds to a large (mx — fJ>x)/T. In that case the first term in ( pip dominates, 

fx ^ exp [(fi x - u a p a )/T] . (32) 

The corrections on account of statistics in ([H]) correspond to an addition (for bosons) or 
subtraction (for fermions) of classical components corresponding to lower temperatures. 
Of a primary concern below will be the classical limit. 

Nonrelativistically, a mean particle kinetic energy for a locally equilibrated system may 
be represented, following (|3~2"D , as consisting of thermal and collective components in the 
form 

(E* n ) = (£f "> v =o + (E c x oU ) = 3T/2 + m x (v 2 )/2. (33) 
4 More elaborate freeze-out conditions are e.g. discussed in |69|| . 



V(±ir +1 exp [n(fix ~ u° Pa )/T} , (31) 
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The collective energy component rises proportionally to particle mass, provided that (v 2 ) is 
the same for different particle species. A relativistic extension of this result may be obtained 
by expanding the relativistic mean particle energy in powers of a collective velocity v, 

(E* n ) = (£f%=o + (E C X M ) « (£jr)v= + ((^>v=o + 2T) (v 2 )/2 + ■■■ (34) 
The mean kinetic energy of particle with mass m in a globally thermalized system, shown 



for reference in Fig. lf| is (E km ) v= Q = 3T + m{K\ — K-i)jKi^ where K u is a modified 
Bessel function of the argument m/T. Relativistically, the thermal agitation increases the 
coefficient in the collective energy component, which is a mass in (|33]), by the kinetic thermal 
energy plus twice the temperature (i.e. at least by 7T/2). That is a large increase for pions, 
and can be significant even for protons or deuterons, at temperatures of interest such as 
displayed in Fig. |[ 

The relativistic enhancement of the collective effect with thermal agitation, suggests the 
possibility of a light-particle distribution that is concave up on account of relativity, when 
the collective motion is present, and a logarithm of the distribution is plotted vs energy. 
On the other hand, a distribution of collective velocities peaked at finite velocity values 
might generate a particle spectrum at low temperatures or for large particle mass, that is 
concave down. The particle momentum distribution is given by 

d 3 N x g x 



dp 3 ~ (2tt) 
and the slope of the distribution is 



drfx (35) 



1 d (d 3 N x \ If E x \ 

lo S \—nr \=m \7> P -— - (7n-v) p , (36) 



T eff dE x °\ dp 3 J T p 
where n is a unit vector in the direction of p and 



(•) p = / drf x (•) / / drf x . (37) 

The slope at p = becomes 
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Fy f = f (^)p=° ~ -jr^( n • v ) 2 )p=o) 



1 ft _l ^ 2 )p=o rn x ((n-v) 2 ) p=0 \ 

TV 2 T ' (38) 



where the approximation leading to the last expression is for low v. Note that the slope 
may be negative, which corresponds to a situation when the spectrum has a maximum at 
a finite p. For a high momentum p, the dominant contribution to the integrals in (|37| ) would 



come from a velocity directed along p to an extent possible given a velocity distribution, 



on account of the strong dependence of (|32|) on v for large p, of a magnitude v m < p/Ex 



assuming a rapidly falling velocity-distribution. The slope may be then estimated as 



T e.f.f T V p 



fm 1 A f 2 



+ — (39) 
TVl + w m TV 2 j v; 

where the second approximation is for Ex ~ p and the last approximation is for small v m . 

If (|38|) is larger than (|39|) then the the spectrum following from ([32]) is necessarily concave 

up within a certain range of energies; if ( |39| ) is larger then the spectrum is concave down (see 

also [pS[|). For example, if collective expansion has a cylindrical symmetry and collective 

speeds are sharply concentrated around a single value v, then the condition for the spectrum 

at 90° with respect to the symmetry axis to be concave down becomes 

2T 

m x > — • (40) 

v 

Clearly, we may have a situation when a heavy-particle spectrum is concave down, while 
a light-particle spectrum is concave up. For a spherically symmetric expansion, the factor 
of 2 on the r.h.s. of (|40| ) gets replaced by 3. Explicit form of a distribution for a cylindrically 



symmetric expansion with a single value of v is (see also ||69|| ) 

d 3 N x ( 1 E X \ T ( 7pvsm0 \ 

~W = M exp {—) Io {—T— ) ' (41) 

where M is a norm, Io - a modified Bessel function, and 6 - an angle relative to the symmetry 
axis. For a spherically symmetric expansion with a single value of v, the distribution in 
momentum is 
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d 3 N x ( 7 E X 

Al exp 



S mh(f9/(f;) 



(42) 



dp 3 V T 

If a chemical potential /ix is uniform throughout a volume V, then the normalization con- 
stant in ([O) and (f£^) is A/" = (7x^exp(/zx/^)/(27r) 3 . For smooth distributions of collective 
velocities, the sufficiently heavy particles (and in principle also possibly the sufficiently en- 
ergetic) probe directly the collective velocity distribution 

d 3 N x g x 



dp 3 



(27T) 

9x 

(27T) 
I 





fx 




exp 




3/2 


nm x / 





/^x — m 



T 



exp 



(2*) 



J dr exp 



m 3 x T 
fjL X -m 



v r 



P 



5 



P 



v r, 



(43) 



where the argument in the exponential in f x was expanded around the maximum, as a func- 
tion of v, and it was assumed that the velocity distribution, nominally defined as 
d 3 V x 



dv 3 



J drdpf x 5(v- v(r)) j J drdp f x 



7jV' / dr exp 



\i x -m 



5 (v — v(r)) 



(44) 



changes slowly over the range in w of a magnitude 7 _2 (T/m) 1 / 2 . 

Besides collective motion, the spectra of different particles may be affected, at low mo- 
menta, by resonance decays and statistical effects. Following (|3T|), corrections for statistics 
in the case of Eqs. fl4"I ) and (^) amount to adding or subtracting from the r.h.s. terms 
corresponding to lower temperatures. In order to relative importance of different 

effects, it will be necessary to understand the dynamics of production. With pions not 
present in the initial state, an important pion observable that may be affected by collective 
expansion, besides spectra, is multiplicity. 



B. Dynamics of Production 



Production of pions in central La + La reactions at 800 MeV/nucleon will be considered. 
At different values of beam energy, pion multiplicities have been measured in La + La 
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reactions ||. At 800 MeV/nucleon, both pion and proton spectra were determined at 
cm. 90° in the very central reactions Using information in ||70|| , these spectra may 

be conveniently normalized for comparisons with calculations. 

Figure [T7| shows the evolution of baryon-density profiles, along and perpendicular to the 
beam axis, in the La + La reaction at b = 0. The system goes through similar stages as 
Au + Au in Fig. [| Figure [18| shows further baryon density at the center of the La + La 



system, net participant collective energy per nucleon calculated as in Eq. (|2lf) , and total 
entropy as functions of time. The collective energy saturates at about 76 MeV/nucleon, out 
of which ~ 56 MeV/nucleon is associated with transverse directions and ~ 20 MeV/nucleon 
with longitudinal. The difference in the energy per degree of freedom for different directions 
is not as large as in the 400 MeV/nucleon Au + Au reaction due to a smaller system size, 
combined with a larger Lorentz contraction and a reduced stopping at the higher beam 
energy. 

One deficiency of the present model is that composites produced during the expansion 
do not participate in pion absorption and production, while in reality they would. To rem- 
edy this, some details relevant to pion production, such as variation with time of particle 
collision-rates and A decay rates in Fig. are presented for a calculation with production 
of composites switched off. However, final pion spectra are given for a calculation with 
composite production switched on, as composites raise the freeze-out temperature, reducing 
to a certain degree the slopes of spectra. The final pion multiplicity does not depend in 
practise on composites as will become apparent later. 



In Fig. 19 it is seen that the rise of A-decay rate is somewhat delayed relative to the rise of 



A-production rate in NN collisions. This is due to a finite A lifetime (see also |73|). The rise 
of the resonance formation rate in Nir interactions follows very closely the rise of A decays, 
on account of a sizeable average Nn cross section. With the exception of the A-production 
rate in NN interactions, that maximizes earlier, all rates maximize at t ~ 13.5 fm/c, when 
the shock waves in lanthanum nuclei reach the vacuum, see Fig. [T7|, and the system as 



a whole becomes excited and compressed. A later decrease of rates in Fig. |17| is associated 
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with expansion. 

As pions (within the model) can only originate from nucleon-resonance decays, and pions 
can only be absorbed into nucleon resonances, the final pion number is equal to the difference 
in the total number of resonance decays and the number of resonance formations in 7rN 
interactions, i.e. the difference in areas under the curves in the bottom panel of Fig. [T9. 



(The N*'s of which the number at 800 MeV/nucleon at any one time never exceeds 4% of the 
total number of resonances, may be ignored in this discussion.) The rather small difference 
between the rates compared to the rates themselves at most times suggests a good kinetic 
equilibrium of pions and deltas with nucleons during the reaction. As, clearly, the number 
of deltas contributing a final pion is identical with the net number of deltas produced in 
nucleon interactions, the final pion number further coincides with the difference in areas 
under the curves in the center panel of Fig. [19]. If one considers a chain starting with a delta 
production in an NN interaction followed by a delta decay with pion emission, followed 
by an absorption of this pion into a delta and so on, then such a chain terminates either 
with a pion in the final state, or with a delta absorbed in an interaction with a nucleon. 
The average number of pion absorptions within any such chain is equal to the ratio of an area 



under the dashed line in the bottom panel in Fig. to an area under the solid line in the 
center panel, i.e. ~1.5. If only chains terminating with a final pion are considered, then the 
number of absorptions is found to exceed 2 in the simulation. One can say that a pion in the 
final state of the system has been, on the average, absorbed more than twice and reemitted 



before reaching the vacuum (see also recent papers |71|j72| ). 

Pion multiplicity measured in central collisions is lower by nearly a factor of 2 than 
anticipated in the fireball model which assumes a complete global equilibration || . A naive 
expectation might be that within the reaction the time is too short to allow for a production 
of the number of pions and deltas expected in a global equilibrium. (The deltas left out 



towards the end of reaction would contribute to a final pion number.) Figure ^ shows the 
pion and delta number divided by total nucleon number, as a function of time in the La + La 
reaction at 800 MeV/nucleon. The number reaches a maximum right after the whole matter 
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has entered the hot compressed region. The ratio of the global-equilibrium number of pions 
and deltas to baryon number at this time, t ~ 14 fm/c, is indicated by the upper of two 
horizontal dashed lines in Fig. |2(]. The ratio is obtained using the equilibrium densities 

dp 1 [dp 1 



(2tt)3 exp(^/T) - 1' PN ^ J (2vr)3 ex P ((E N - p N )/T) + 1 



f dp f dm 1 
PA = 16 JJ2^J 2^ A exp((£ A -M/T) + l' (45) 

where temperature and chemical potential are adjusted to reproduce the representative 
baryon density, p = 2.25p at a respective time (cf. Fig. |i~7|) , and the energy per nucleon 
from simulation, excluding potential energy, of 1133 MeV/nucleon. (Isospin asymmetry 
is ignored for simplicity.) The equilibrium value, indicated by the upper dashed line in 
Fig. |20], is not reached, what appears to support an expectation that the time is too short 
to produce an adequate number of pions and deltas. In fact, as the system expands and 
baryon density drops, the global equilibrium value increases more and likewise the apparent 
deficit of pions and deltas. Towards the freeze-out, however, the number of pions and deltas 
behaves in an unexpected manner. Instead of straining towards the equilibrium value and 
settling once interactions stop, this number declines. The decline is seen in many other 
microscopic calculations of collisions, e.g. in |66| , |73H 75|] . (Incidentally, if a ratio of pion and 
delta density to baryon density is investigated locally, then the decline becomes even more 
pronounced; the total number represents the system as a whole with different parts going 
through different stages.) The rise and decline in the number of pions and deltas could be 
understood if somehow the number of pions and deltas were first below and later above the 
equilibrium number. The clue to an understanding of this situation is in collective motion 
ignored till now. 

Collective motion within the system takes away energy available locally for particle 
production. Right after the whole matter in the La + La reaction has entered the hot 



region, at t ~ 14 fm/c, the collective energy within this system, as seen in Fig. [18|, 
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is ~ 30 MeV/nucleon. If this energy is subtracted from the available energy and Eqs. fl45|) 
are used as before, then the equilibrium value for the ratio of pions and deltas to baryons 
is obtained which is indicated by the lower of two horizontal dashed lines in Fig. It is 
apparent in this figure that the system approximately reaches an equilibrium with regard to 
the pion and delta number given the particular value of collective energy. As the sys- 



tem expands, though, more and more energy turns into collective. In Fig. |21] it is 
seen that if expansion were isentropic then the number of pions and deltas should de- 
crease. (At low temperature, in an isentropic expansion the pion-to-baryon ratio behaves as 
Ptt/p oc exp(— m n pl^ 3 /Tip 2 / 3 ).) In Fig. |2l], the state of the system at t — 14 fm/c is crudely 
(as in fact different parts of the system go through different stages) represented with crosses. 
To assess the situation further, the equilibrium number of pions and deltas is calculated at 
different times by evaluating the equilibrium density of particles throughout the space in 
local frames, given energy per baryon and baryon density in such a frame, and integrat- 



ing the particle density over whole space. Results are represented with dots in Fig. |20 . 
The obtained local-equilibrium number of pions and deltas is seen to increase when the 
mass of hot matter increases, and to decrease rapidly afterwards. The actual number of the 
particles closely follows the equilibrium number during compression. The local-equilibrium 
number crosses the actual number of pions and deltas at about the maximum of the latter 
number (as should be expected). While the number of pions and deltas decreases during 
the expansion, it freezes out, nonetheless, at a relatively high value. This is because a de- 
crease of nucleon density in expansion effectively shuts off the processes N + N ^ N + A, 



see Fig. [19|. The processes of pion absorption, with a larger cross section, and the A decays 
continue somewhat longer, preserving the thermal but not the chemical equilibrium of the 
pionic sector with nucleons. Portion of the chemical equilibrium path, with nonadiabaticity 
taken into account, is sketched with dashed lines in Fig. |21]. At a freeze-out, at a density 
p ~ 0.5 po, the ratio of pions and deltas to baryons is larger by ~ (1.4 — 1.8) than in the 
state of chemical equilibrium. 

It has been conjectured in the past || that pion multiplicity in the central collisions is 
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low, because part of the energy in the center of mass is used up in compressing matter, 
and, correspondingly, is not available for production. Present results show that the role 
of energy not available for production within a heavy system is played by the collective 
energy In fact, potential energy at the time when pion and delta number maximizes in 
the La + La reaction simulation, differs from the potential energy in a ground state by 
a mere 3 MeV/nucleon. (Momentum-dependent interactions might, though, have an effect 
on pion number.) With regard to another expectation, that the time in a reaction is too 
short for producing an equilibrium number of pions and deltas, it is seen that the number 
of pions and deltas actually freezes out not below, but above an equilibrium value. 



C. Features of Single-Particle Observables 



The left panel in Fig. |22] shows calculated momentum distributions of different parti- 
cles (filled symbols) within a cm. polar-angle range of 60-120° in the 800 MeV/nucleon 
La + La reaction at b = 1 fm, together with momentum distributions of protons and neg- 



ative pions in the very central reactions deduced from experimental results of Ref. fTO 



(open symbols). The cross section for central reactions investigated in Ref. [ 70fl , needed for 
the normalization of distributions was estimated at 86 mb, using the integrated 7r~ yield 
stated in ]7(| and tc multiplicities from Ref. ||. The experimental tc~ distribution agrees 
very well with the calculated distribution at all energies, and the proton distribution - up 
to an energy of ~ 450 MeV in the c.m.s. The pion distributions may be well approximated 
with straight lines (in the logarithmic plot), while baryon distributions are concave down. 
The slopes significantly decrease with increasing particle mass. If the calculated slopes were 
to be interpreted in terms of temperature, quite disparate values would be obtained, of 70, 
~ 110, ~ 145, and ~ 175 MeV, for pions, protons, deuterons, and helions, respectively. 
For the sake of understanding the momentum distributions from simulations, the left 



panel in Fig. 22, shows additionally baryon and neutral-pion distributions from an instan- 



taneous freeze-out model (solid lines). In the case of an expansion with a single value of 
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velocity, the asymptotic collective energy in Fig. [18] would correspond to v ~ 0.41c. In the 



case of a uniform velocity distribution, adopted within the model for the display in Fig. |22|, 
the collective energy corresponds to a cut-off velocity- value of v c ~ 0.51c. The mean proton 
energy at a freeze-out of 196 MeV from simulation, corresponds then to a local temperature 
of T « 55 MeV. The baryon momentum distributions within the model take the form 

^ - g (f ) If- 2 - (-^) - (¥9 /(¥) . w 

cf. (|2]). For pions the distribution is (see also |76| , |77| ) 
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where the second and third term inside braces account, respectively, for the effects of Bose 
statistics in lowest order, and for the decays of A resonances left-out at a freeze-out. The min- 
imum and maximum momenta of a A whose decay might yield a pion with a momentum p^, 
are (pa)i,2 = ( m / m V)\^P^ T EnPt\- The subscript '*' indicates pion momentum and en- 
ergy in a A-frame. Pions and deltas are assumed to be in a thermal but not in a chemical 
equilibrium with the rest of the system; the chosen value of a freeze-out density is p = 0.6po- 
It is quite apparent in Fig. ^ that the instantaneous freeze-out model (with parameters 
taken from dynamics) describes, in a semi-quantitative manner, the shapes of spectra from 
simulations. The spectra of the lightest particles, pions and protons, with largest thermal 
dispersion of velocities, are, in fact, equally well described within a model with a single 
collective expansion velocity. Proton mass (and cluster masses) satisfies m > 3T/v and the 
spectrum is concave down. For pion mass the opposite inequality is true. Long- and short- 



dashed lines in Fig. |22| show, respectively, contributions from free pions and from A decays 
within the instantaneous freeze-out model. The free-pion contribution is concave up, but in 
the overall distribution this concavity is partially filled up with pions from decays. Despite 
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the spread of the A resonance in mass, a major resonance contribution comes here (and at 
other temperatures of interest) from masses in the vicinity of resonance peak; the ratio of 
resonances to pions is pa/ Pit oc pn exp(— (r/iA — — m^) /T) /T 3 ' 2 . 

Of a general interest is the degree of Bose enhancement or the proximity to a conden- 
sation for low-momentum pions. The excessive number of neutral pions at a freeze-out 
corresponds to a finite chemical potential /i^o ~ T\og{N n o /N^q) ~ 25 MeV. For nega- 
tive pions, the isospin asymmetry within the system corresponds to an additional enhance- 
ment of the potential for negative pions, and a depletion for positive pions, fi n - — fi n o = 
p, n o — fi w + ~ T\og(N/Z) ~ 20 MeV. Overall, the potential at a freeze-out for negative 
pions may be as large as ~ 45 MeV, corresponding to a low-momentum enhancement of 
exp((/i 7r - — m 7r -)/T) ~ 18%. Possible enhancement for neutral pions is 13%. Beyond that, 
the Coulomb field in the final state enhances distribution of negative pions at low momentum 
and depletes that of positive pions and other particles. Lower limit on the magnitude of the 
effect on single-charged low-momentum particles follows from the formula obtained under 
an assumption of a dominating collective motion, 0.7 MeV-fm p l l 3 Z / A l l 3 {E km ) ~ 5%, where 
Z and A are total charge and mass numbers, respectively, (E km ) is mean kinetic energy of 
particle species in question, and p is freeze-out density. Inspection of spectra reveals that 
the actual effect of the Coulomb field in the La + La reaction is of the order of 30%. 

Collective motion is generally characterized by a correlation between particle spatial 
positions and their momenta. In order to demonstrate explicitly that strong features in 
the calculated distributions, particularly slope-differences between different particle species, 
are associated with the collective expansion, calculations are further repeated with an inter- 
change of particle positions during evolution. Interchanges are made within same species, 
when collisions are still frequent. (Particle momenta are not altered.) The procedure leaves 
density in configuration space and momentum distributions intact. Correspondingly, kinetic 
and potential energies stay the same. However, correlations between particle momenta and 
positions get destroyed. The final momentum distributions from calculations with the po- 



sition interchange are shown in the right panel of Fig. Straight parallel lines are drawn 
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to guide the eye. The distributions from the calculations turn out largely exponential with 
nearly the same slopes (some collective motion still develops due to the collisions, cf . Fig. |18|) . 
There is enough time within a reaction for a complete thermalization. 



Figure |23] shows measured pion multiplicities in central symmetric reactions as a func- 
tion of beam energy [[|, together with multiplicities from standard calculations of central 
La + La reactions (solid line) and from calculations with positions interchanged (dashed 
line). The standard calculations, with the collective motion developed during evolution, de- 
scribe reasonably the data. The calculations with the position interchange and the collective 
motion extinguished, on the other hand, lead to multiplicities larger by nearly a factor of 2 
(essentially such as expected within the fireball model), in a dramatic disagreement with 
data. 

Despite of the fact that the measured pion multiplicities point strongly towards the col- 
lective motion, in order to better assess this motion, it would be clearly highly desirable to 
have normalized pion and proton momentum distributions from central collisions of heavy 
nuclei, extending over a wider kinetic-energy range than in Fig. [22|, possibly together with 
cluster distributions. Heavy systems are preferred over light systems [f|]. In peripheral col- 
lisions of light systems, differences between pion and proton spectra may exceed differences 
arising from collective-energy considerations, on account of energy conservation in elemen- 
tary collisions, with kinetic energy depleted due to the production of pions. In principle, 
with a decrease of impact parameters, the relative importance of conservation in individual 
collisions would decrease, while the importance of collective motion should increase. To test 
an interpretation in terms of collective energy, features of the cluster distributions could be 
examined | 78[| . In lightest systems, though, such an interpretation would be always ham- 



pered by the lack of a complete stopping even in most central collisions, at beam energies 
at which pions are abundant. 

The actual distribution of collective velocities for baryons in the vicinity of last collisions 
from the simulation of the 800 MeV/nucleon La + La reaction is shown in Fig. |2I|. This 
distribution is not peaked at any large value of velocity (although v 2 (d 3 V/dv 3 ) is peaked, 
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however broadly, in the vicinity of v ~ 0.4c). Tests indicate that the high- velocity tail is 
associated with low temperatures. High temperatures, on the other hand, are particularly 
associated with the low-velocity portion of the distribution. Parametrization of particle 
spectra, such as from experiment, in terms of a single freeze-out temperature for any particle, 
will be always approximate. Also an approximation can generally arise from the use of same 
temperature and same velocity for different particles, as e.g. pions may be emitted earlier 
when temperatures are slightly higher and the collective energy lower, than when light 
baryons are emitted, and differences generally develop in the emission of different baryons, 
cf. Sec. [VI|. In any case, though, the use of a uniform distribution in collective velocity in 
parametrizing particle distributions would be more sensible than the use of a single velocity 
value. One possibility within an experiment is the use of a velocity distribution deduced 
from heavier fragments, of which distributions should approach that of the collective velocity, 



with increasing fragment mass, cf. Eq. (f43|). The use of a single velocity value and of 



a temperature treated as adjustable parameters in describing measured spectra of different 
particles, on the other hand, is expected to give absurd results. Reasonable temperatures 
would lead to excessively sharply-peaked distributions of heavy fragments. To make the fits 
acceptable, temperature parameters would then need to be raised making it appear that 
heavier fragments came from hotter regions than light fragments. 

D. System Measured by Brockmann et al. 

Before moving on to other issues, let me discuss data ||, that first stirred up interest 
in shapes of transverse spectra of pions from central reactions, from the Ar + KC1 system 
at 1.8 GeV/nucleon, and further aspects of transport models. In [[J experimenters attempted 
to fit their data with a single exponential and concluded that this was not possible. Moreover, 



they compared results to those from a cascade model |f75| . Due to the low mass of the 
system and high beam energy, and corresponding lack of stopping, the spectra from transport 
simulations are quite sensitive to the assumptions on elementary collisions; the experimenters 
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chose a cross section for N + iV — > N + A processes that was exceedingly forward-peaked. 
Cascade results were not laid over data, but rather fitted with an exponential. The conclusion 
was that the cascade model cannot explain the /jzg/i-momentum portion of a pion spectrum. 

The data from || are shown in two panels of Fig. |25] (open circles). Insets show an asym- 
metry coefficient a obtained by fitting the double-differential distribution 

d 2 AL 



dE d cos 6 



oc (I + a(E) cos 2 6). (48) 



The right panel shows further results from the cascade model [[F3|] following assump- 
tions in || (squares). It is seen that, indeed, the high-momentum tail of the measured pion 
distribution is underpredicted. Further, the asymmetry coefficient from the model is far too 
high at high pion kinetic energies. Discrepancy in this coefficient reveals part of the problem 
with the high-momentum tail. The left panel of Fig. [25| displays, in addition to data ||, 
results from the present model |35| with an angular dependence of the delta-production 
cross-section that was fitted to elementary data in Ref. |0J] (filled circles). The model de- 
scribes well the high-momentum tail of the pion spectrum, and it describes much better 
the asymmetry coefficient than the cascade model. This is associated both with the shape 
of elementary delta-production cross-section in the model and with respecting of the time- 
reversal invariance in absorption cross-sections (cf. Sec. U high-energy pions are excessively 
absorbed in the original cascade model compared to the low-energy pions). For illustrative 
purposes (the involved assumption is too extreme), calculations were repeated within a ver- 
sion of the cascade model where A-production cross-section was isotropic; the results are 
shown in the right panel of Fig. |25] (diamonds). The description of the high-momentum tail 
and asymmetry coefficient is improved over the first cascade calculation. 

Overlaying the calculations and data reveals a puzzle concerning the /owMnomentum 
part of pion spectrum. Already the cascade model underpredicts data |J in this region 
and correcting the pion absorption cross-section, cf. Eq. (||) and associated text, increases 
discrepancy. Measured spectrum appears to have a low-momentum peak as though super- 
imposed onto the spectrum such as from the calculation, of a similar height as the latter, see 
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left panel in Fig. |25], apparently pointing towards physics outside of the transport model, f] 

E. Shadowing 

As impact parameter is increased in symmetric collisions at high beam energies, regions 
with spectator matter develop and grow in size. In strongly asymmetric systems, such as 



investigated in [p0|| , the spectator matter is always present. The matter can shadow pions 
emitted from the hot participant region. 

In Fig. ^ with mean momentum-components within the reaction plane of particles emit- 
ted from 1 GeV/nucleon Au + Au reactions, it is seen that, at low impact parameters 
~ 3 fm, pions moving forward and backward in the c.m.s. have momentum components 
of same sign as baryons. As the pion angular-distribution is forward-backward peaked at 
b = 0, cf. Fig. |7|, the average deflection such as for baryons, within the reaction plane, must 
be associated with the weak discontinuity. The pion momentum components per mass are 



considerably lower than the baryon components, see Fig. [L3[ and pion velocity components 
even more so, by a factor of ~ 5 in the forward and backward rapidity regions, than the 
baryon components. At higher impact parameters (~ 7 fm) in the reaction, pions have 
momentum components of an opposite sign than baryons, cf. Fig. [13|. These components in- 
crease more than twofold if analysis is limited to pions emitted early, within the first 16 fm/c 
of the reaction. (The latter pions constitute about 1/3 of all.) 

When pion emission is examined at a 90° cm. polar-angle, it is found essentially isotropic 
in the azimuthal angle at low impact parameters such as b ~ 3 fm. However, at higher 



5 Elsewhere, though, results well exceeding data in the low- momentum region were obtained [79], 
supposedly on the account of the nuclear mean field and Pauli principle. Here, little effect of these 
factors on spectra is found (as is apparent in Fig. |25| ) in accordance with the intuition, given the 
small system size, the high beam energy, and the momentum region. No attempts are made here 
to interpret unnormalized data in literature. 
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impact parameters, such as b ~ 7 fm, this emission is found peaked at 90° with respect to 
the reaction plane. The peaking is enhanced when analysis is limited to pions with high 
cm. momenta, p > 200 MeV/c, see Fig. ^6] (G = 0°). The latter peaking was observed 
experimentally in events corresponding to intermediate impact parameters |j5I^52| . 



Figure displays contour plots of baryon density (solid lines) within the reaction plane, 
in the 1 GeV/nucleon Au + Au reaction at b = 7 fm, at a time t — 15 fm/c. At the impact 
parameter equal to the half of the maximum possible value, the spectator regions are well 
developed. The spectator matter begins to detach from participant matter at the time at 
which the densities are displayed. Shocks have started out in the particular reaction at 
an angle of inclination relative to the beam axis of ~ 65°, cf. Eq. fl25|), and have rotated 



down to lower angles ~ 55°, compare Fig. 11. The last of what has become the participant 



matter has dived into the shocked region at a time t ~ 10 fm/c. 



Figure |27] displays further contour plots of pion and delta density within the reaction 
plane. It is apparent in that figure that the spectator matter limits pion emission. An analy- 
sis of the emission in different directions, within the simulation, shows that the transmission 
coefficient for pions through the center of spectator matter does not exceed 5% during the 
first 16 fm/c of the evolution. With not very strong sideward-flow effects exhibited by pions 
in the absence or near absence of spectator matter at low impact parameters (at this and at 
other beam energies), the evidently strong shadowing at the higher impact parameters must 
generate sizeable average momenta per mass opposite, due to geometry, to the momenta of 
nucleons. 

Early proton emission is affected by shadowing, too. In the proton case, though, the shad- 
owing does not lead to such dramatic effects as in the pion case, because baryons exhibit 
generally stronger dynamic effects and because baryon emission is overall somewhat delayed 
relative to pion emission. Thus, the mean momenta within the reaction plane of protons 
decoupling from the 1 GeV/nucleon b = 7 fm Au + Au, within the first 16 fm/c, are nearly 
zero. These protons constitute a smaller fraction of particles emitted at all times than pions 
(about 10% as compared to 1/3 for pions; pions get ahead of the bulk of the matter on 
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account of larger thermal velocities). 

While the effect of shadowing on the mean momenta is opposite to that of dynamics, 
both type of effects act to suppress the emission within the reaction plane at midrapidity 



as compared to the emission out of the reaction plane. Figure ^6] with calculational results 
illustrates the fact that the issue as to which of the above two effects dominates in gen- 
erating an azimuthal pattern at midrapidity, may be resolved at a qualitative level within 
an experiment. Thus, one needs to examine an emission pattern within a plane p z> = 
where axis p z ' points at an angle G measured counterclockwise within the p x -p z plane, from 
the beam. When shadowing dominates, an azimuthal anisotropy should persist or increase 
compared to that at = (within a plane perpendicular to the beam), when axis p z> is ro- 
tated to negative 0-values, making the plane of analysis more aligned with spectator pieces 
when latter are close to the participant matter, at a time when many particles are emitted, 
cf. Figs. ^7] and §. On the other hand, an anisotropy solely due to shadowing should decrease 
or disappear when axis p z> is rotated to positive B-angles (case of pions in Fig. ^6|). When 
anisotropy is primarily associated with a stronger collective motion in the directions parallel 
to shock-front surfaces, in particular in the direction out-of-the- reaction-plane, in compari- 
son to the direction perpendicular to shock surfaces, then the anisotropy within the plane 
p z> = should increase when axis p z> is rotated to positive G-angles (case of protons in 



Fig. and baryons in Refs. [ 13|,p7 |). 



With regard to the transverse-momentum dependence of azimuthal anisotropies [pl| , p2 



any anisotropy must disappear whenp- 1 — > (but see ||71||), provided that particle momentum 
distribution, denoted below as 
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is smooth at p L = 0. Expansion of the distribution up to the second order in transverse 
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(49) 



where is azimuthal angle relative to the reaction plane and x and y denote directions 
perpendicular to the beam axis in- and out-of-the-reaction-plane, respectively. The coeffi- 
cient in front of cos vanishes at midrapidity in a symmetric system. The coefficient in front 
of cos 20 is proportional to momentum squared. (These features of the coefficients would 
be retained if third axis were rotated away from the beam.) The azimuthal asymmetry 
at midrapidity may be e.g. quantified with a ratio of out-of-plane to in-plane momentum 
distribution minus 1. This quantity is proportional to the coefficient in front of cos 20 in 
Eq. (^) and, correspondingly, vanishes quadratically with momentum as p — > 0. (That 
is also true for other measures of asymmetry such as (cos 20).) In the case of pions the 
asymmetry would grow with an increase of pion velocity at midrapidity, as pions would 
move out to the sides in a reaction and encounter spectator matter in the vicinity of the 
reaction plane. At lower impact parameters, when layers of spectator matter are thin and 
semi-transparent to pions, some dependence of the asymmetry on pion momentum can be 
expected, related to the energy-dependence of pion-nucleon cross-sections. At high impact 
parameters in heavy systems, when layers of spectator matter are thick and intransparent 
to pions even for low pion-nucleon cross-sections, the asymmetry should saturate once pion 
velocities get close to c. When this paper was under completion, an extended study of 
pion shadowing by Li has appeared ||83|| , with calculational results exhibiting these features. 
Data |y],|82| are suggestive of the rise of asymmetry at high momenta, but corresponding 
errors are large. Other recent theoretic references on pion shadowing in symmetric systems 
are U7HJSTO 



VI. BARYON OBSERVABLES 
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A. Entropy and Fragment Yields 



Experimentally yields of different fragments are used to determine entropy produced in 
reactions. Specifically at low densities when effects of statistics and interactions between 
fragments can be ignored, thermodynamic entropy per baryon can be expressed as 

A ^ A \2 TJ~2AATAT 1 ' 

where Mx is multiplicity of fragment X, Bx (> 0) is fragment binding energy, and fix = 



Zx/J>p + Nx/^n is chemical potential. The approximation in (pOj) is for small binding energies 
per nucleon in comparison to nucleon chemical potentials fi n and \i v . The ratios of nucleon 
chemical potentials to temperature can be obtained from fragment yield ratios, thus e.g. 



iu n /T = \og(g p m 3 p / 2 M d /g d m 3 d /2 M p ) = \og(M d /3V2M p ), and fi p /T = \og(^2/3M h /M d ). High 
entropy values correspond to large negative values of Hn/T and to a rapid drop of yield 
with fragment mass. In data analysis |2y|l^|, theoretic models are used, such as QSM ||86|| , 



to extend the relation (|50| ) between entropy and fragment yields to situations when e.g. 
fragment phase-space overlap becomes of a concern. 

Entropy per nucleon determined from measured yields of light fragments at large cm. 
angles in the Au + Au reaction at 400 MeV/nucleon |2{]|] , is shown in the left panel of Fig. ^ 
as a function of reduced participant-proton multiplicity. The right panel of Fig. |28| shows 
entropy per nucleon associated with fragments emitted at large angles |35| (circles), deter- 
mined from calculated fragment phase-space distributions, as a function of b. An agreement 
is found with regard to the value of entropy produced in the very central reactions and with 
regard to the general behavior of entropy with the reaction centrality at the particular beam 



energy; see further Ref. [pq] . 



In order to understand the production of entropy in collisions it is convenient to consider 
first the case of b = 0. In that case, at high energies, most of the entropy is produced in 
shock fronts, cf. Figs. [5] and f|. Further production takes place during expansion. The hy- 
drodynamic Navier-Stokes equations yield for the rate of change of entropy per baryon with 
time |^0| 
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^UJ +V - V UJ = ^ dlV(KVT) + 2^(^ + ^-3^^J ' (51) 
where k is heat capacity coefficient and rj is viscosity coefficient. In the classical limit, 



cf. e.g. pBfl , the coefficients are approximately given by k ~ (75^/64<jNN)yT/m and 
rj ps (5y/Tr/16)y/mT/aNN- However limited, an insight into the entropy production dur- 
ing the expansion may be obtained by substituting into fl5"T| ) results for velocity, density, 
and temperature from the simplified consideration in Sec. [Ill B| , involving assumptions of 
cylindrical symmetry, self-similar expansion, and ideal-gas equation of state. An examina- 
tion of the r.h.s. of Eq. ( pi]) shows that the heat conduction can cause a drop of local entropy 
with time when temperatures on the average decrease with moving away from a given point, 
whereas the viscosity always leads to a rise in the local entropy per particle. Substitution 
of the results from Sec. |I[B] yields a net drop of entropy at the center of the system with 
time, quantitatively consistent with what is observed in Fig. |3|; temperatures decrease in 
the surroundings due to the expansion into the vacuum. The substitution further yields 
a divergence of entropy production-rate towards the edges of density distribution; in the net 
entropy production-rate the viscous contribution is integrable while the conduction contri- 
bution is nonintegrable. If a cut-off density is adopted for the validity of hydrodynamics, 
then the net increase of entropy per nucleon on the account of heat conduction in expansion, 
may be estimated with the formula 

m 3( T - i)V /V(o,o)V /2 

\aJ a NN p (o,o) r{o) \ pi j ( ' 

cf. Eq. ([T7|) . The increase due to viscosity turns out to be given by a similar formula but 
with a numerical coefficient in front smaller by one order of magnitude. Rise of the entropy 
production-rate towards the edges of density distribution is quite consistent with what is 
observed in Fig. The results indicate that most entropy in expansion at b = is produced 
within an interface between the hot matter and the vacuum and that the production is 
associated with the equalization of temperature. Quantitatively, Eq. (|5*2"| ) gives an increase 
of the entropy during an expansion in the 400 MeV/nucleon reaction by 5(S/A) ~ 1, while 
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the simulation yields S(S/A) ~ 0.6. 

Some limitation of the above consideration becomes apparent when a self- consistent so- 
lution of the viscous hydrodynamic equations is confronted with simulations. While finding 



that the consistency in the solution moderates an increase of the entropy, Kapusta [57 
obtained a near independence of the entropy generated in expansion on the mass of an ex- 
panding nuclear system, S(S/A) ~ 0.3 for p\ ~ 2po- By contrast, in the simulations the 
generated entropy increases when A decreases, actually in nearly such a manner as predicted 
by Eq. (j52|) , 5(S/A) oc A~ 1 ^ . The increase appears, in fact, consistent with the central data 
for different systems p9t^0 



Following (pO"D or more detailed considerations, generation of entropy during expansion 
in the considered 6 = Au + Au reaction, increases local kinetic energy towards freeze-out 
by ~ 50%, limiting the rise of collective energy. In principle, the dissipation could have been 
so strong that no collective energy would have been generated at all. Then, correspondingly, 
cluster yields would have been lower also. The correlation between cluster yields and collec- 
tive motion is illustrated in Fig. ED] which shows the beam-energy dependence of yield ratios 



from standard calculations and from calculations with the collective motion extinguished by 
particle-position interchange, together with data [ Z0| . The beam-energy dependence of the 



ratios of measured A = 3 cluster yields to proton yield is not well reproduced within the 
standard calculations; the quality of the description is about the same as within the QSM 
model |2(J. Unambiguously, though, the data favor the results from the standard calcula- 



tions with a significant collective motion developed, over the results from calculations with 
the motion extinguished. 

Dissipation heating the edges of matter distribution and cooling off the center of a col- 
liding system at b = as in Fig. |3|, beyond what actually can be seen there, together with 
some emission before equilibration, lead to quite a spread of entropy values per baryon 
in the vicinity of particle emission points, cf. Fig. |30[ As, generally, more clusters are 
emitted when entropy is low than when it is high, the average entropy in the vicinity of clus- 
ter emission-points may be expected lower than in the vicinity of nucleon emission-points. 
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Specifically, if the distribution such as in Fig. [30 were approximated in a gaussian form 



dV/d(S/A) ~ exp — {S/A — S/ A) 2 / 2 A 2 {S/A) , and a complete equilibrium were assumed 
with Mx/Mn ~ exp [— {Ax — 1)(S/A)] for a given entropy, then the average entropy in the 
vicinity of emission points of baryons with different masses would be expected to decrease 
with a decrement per unit mass (S/A)a+i — (S/A)a £ — A{S/A)/2. In the 400 MeV/nucleon 
Au + Au reaction the dispersion is A(S/A) ~ 0.40, and the expected decrement would then 
be (S/A)a+i — (S/A) a £ —0.20. In the simulation, the entropy in the vicinity of last collision 
points is found to decrease by about 0.14 per mass unit for light baryons. With an increase 
of the beam energy, the entropy generated in expansion increases, just as does the entropy 
generated within shock fronts. Likewise, the dispersion of entropy values increases, and the 
mean entropy in the vicinity of emission points is found to decrease faster with baryon mass. 
E.g. in the 1 GeV/nucleon b = Au + Au reaction-simulation the decrease is about 0.20 
per mass unit. The entropy differences diminish, on the other hand, when the beam energy 
is lowered (though, they become then increasingly more difficult to assess reliably within 
the calculation). 

In the analysis |B| of intermediate-mass fragment yields from central Au + Au reactions, 
in particular at 400 MeV/nucleon, a gradual decrease of entropy with charge was observed, 
however slower than in the present calculation. For heavy fragments the decrease might be 
moderated by size effects. In the analysis |19| the entropy values at the low-Z end were 
further found lower than in [20[, by 0.5 at 400 MeV/nucleon, which could become clarified 
once masses get resolved and angular coverage extended within that experiment. 

At a finite b, generally, more variation of entropy values throughout the system is ex- 
pected than at b = 0, as spectator regions are characterized by S/A ~ (1— 2). With regard 
to the participant region, a shock at a finite b would lead to a lower entropy than a shock 
at b = 0, see Eq. ( p5|) and Fig. |j. If one assumed that about the same amount of entropy 
were produced in an expansion at a finite b, as at b = 0, then one would expect, on the basis 
of the RH equation, a change of entropy with b to follow 
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(I) 6 - (£) 6 o = j (Pi (0, J , - | (Pi (0, E beam ) ,E), (53) 

where entropy on the r.h.s. is written as a function of local baryon density and of net local 
energy per nucleon. The density p\ is from Eq. (|26|) and energies E and E' are related, 
respectively, to E beam and E' beam , cf. Eq. ([25]), according to £ = m N - 16 MeV+£ 6eam /2(l + 



1 + E beam /2mN)- The result for the entropy as a function of b from (|5^) is shown in the 
right panel of Fig. [28| with a dashed line. The decrease is clearly in contradiction with 



the simulation and with data. Note that the experimental entropy |20| is determined from 
fragments emitted into wide angles, i.e. stemming from the participant region. The consid- 
eration, though, ignored till now the weak discontinuity. As the latter spreads, the viscosity 
slows down the nuclear matter, generating entropy. If one assumed that kinetic energy were 
completely thermalized for a given density, then one would expect for the matter stemming 
from the center, instead of (|53|) , 

(j) 6 - (~) fe=o = | (Pi (0, O ,E) - | (pi (0, E beam ) , E) . (54) 

This now leads to an increase of the entropy with an increase of b as the thermalization 
for a given energy occurs at a lower density. The increase is quadratic in b/b max for low 
values of the reduced impact parameter, cf. Q26D, in a qualitative agreement with the data. 
The difference in entropy between high and low b increases with an increasing beam energy, 



as seems to be the case with data f20fl . Quantitatively, though, the entropy from (|54|) is 
still too low compared to the large-multiplicity data displayed in Fig. ^8| or the simula- 
tion at large b. (Although at a maximum compression the entropy gets, indeed, close to 
(S/A) (pi (0, E' beam ) ,E).) The result from (|54"D is shown by a solid line in Fig. ^ and only 
the difference between the solid and dashed lines can be attributed to the viscosity. The 
additional entropy at large b is associated with the rise in the entropy generated during 
expansion, when the size of participant region diminishes. For b close to maximum the 



participant mass is |88[ A part ~ 3y/2A(l — b/b ma x) 2 , where A (= 197) is the mass of one 
nucleus. With the scaling in the simulations such as in (^), the generated entropy behaves 
then as S(S/a) oc (1 — b/b max ) 2 ^ 3 for large b. 
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Summing up the discussion on entropy production, most entropy at b = is due to 
shock waves and some is generated in expansion. As b increases, the shock contribution to 
the entropy in the participant region decreases, and the contribution of the spreading weak 
discontinuity rises from 0. Given that the density at the center gets close to pi from Eq. ( |26| ) 
and entropy close to (S/A) (pi (0, E' beam ) ,E), it is apparent that properties of matter in the 
compressed state at the center can be regulated by selecting impact parameters and changing 
the beam energy. The entropy per nucleon generated in expansion increases with b. Even 
at b = and a high energy, the entropy in the vicinity of particle emission points is spread 
out around the mean entropy in a reaction, because of the dissipation. 



B. Collective Energy 



Very recently, values of mean baryon energies |T5|--|i~8f and shapes of spectra p~T| , [T8|j2T|1 
have been used to extract the magnitude of collective energy in collisions. (Procedures 
followed, in particular, Eq. (|3^) and nonrelativistic versions of Eqs. ( f4~2"|) and (|46|), with 
Coulomb corrections; see also [0].) It is then worthwhile to discuss the experimental results 
on the energy and compare them to calculations. 

The collective-energy values of about 18, 31, and 51 MeV/nucleon, deduced by the 
FOPI collaboration using intermediate mass fragments (IMF) |l7j from central Au + Au 
reactions at beam energies of 150, 250, and 400 MeV/nucleon, come close to the values of 
collective energy at freeze-out in present simulations, of 21, 35, and 55 MeV/nucleon at the 
respective beam energies, see Fig. [|. Revision |18| of the result for 150 MeV/nucleon, with 
a lower bound on collective energy put at 10 MeV/nucleon, stemmed from a concern that 
studied events were not exactly central, and the value of collective energy was inflated due 
to a persistance of longitudinal motion and a particular angular range of the detector setup. 
It is argued below that latter bound is likely too cautious. 

As the separation of kinetic energy into collective and thermal or excitation components 
occurs at freeze-out, this separation should not depend very strongly on impact parameter 
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(cf. discussions in Sec. |T| and Ref. [^2|; excitation energy is not expected to vary rapidly). 



The collective energy per nucleon for light baryons from simulations, as a function of im- 



pact parameter, is illustrated in Fig. [Jl]. Between 6 = and b = 12 fm the energy rises 
only by 20%. For heavier particles, and if statistical decays of spectator fragments were taken 
into account, the rise would be generally faster. Nonetheless, in the {p ± -y) distribution pre- 



sented in |L8[ no traces of any spectator decay can be seen. An average impact parameter is 
{b) < 3 fm, and particle charge is 3 < Z < 8. Relative to b = 0, for such charges and impact 
parameters, still only a small change of collective energy would be expected. What definitely 
changes with the impact parameter is the division of collective energy between the longitu- 
dinal and transverse degrees of freedom. If simulations adequately described the situation 
in reactions, then, at representative impact parameters, the energy per transverse degree of 



freedom would be close to the longitudinal energy, and the angular cuts in [0,0 would not 
affect the deduced value. (Note, that the equality of energies per degree of freedom does 
not contradict azimuthal anisotropies ||18|| , cf. Sec. LV B[ ) Given an unfavorable situation for 



the determination of collective energy, with an anisotropy similar to that in the simulation 
at b = 0, but with more collective energy associated with the longitudinal degree of freedom 
than with one transverse, one could still put a more stringent lower bound on the collective 



energy than in ||18|| , at 14 MeV/nucleon, taking into account angular cuts and a behavior of 



particle energies with Z, and assuming smooth momentum-distributions. Analogous lower 



bounds, about 25% below the initially deduced values [0, at 27 and 40 MeV/nucleon, re- 
spectively, could be put on the collective energy in the central Au + Au collisions at the 
two other beam energies of 250 and 400 MeV/nucleon. 



At b ~ 0, or when examining transverse directions one might expect a gradual 
decrease of collective energy with mass. Thus, in Fig. |3] it can be seen that regions with 
a high entropy, from which few heavy particles would originate, freeze-out early and are 
associated with relatively high collective velocities. Regions with low entropy, on the other 
hand, would freeze-out late and would be associated with lower than average collective 
velocities. Spectra of IMF, corresponding to wide angles in the cm., from Au + Au reactions 
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at 100 MeV/nucleon and b < 4 fm, were analysed in [[21]] , and collective energy per nucleon 
gradually decreasing with fragment mass was deduced. The mean collective energy at the 
low mass end was ~ 9 MeV/nucleon, which should be compared to (cf. |[22|| ) (3/2)E t 
11.5 MeV/nucleon in the simulations of that system at the respective impact parameters. 

VII. DEPENDENCE ON INTERACTIONS 
A. Collective Energy 



Judging from Fig. glj little sensitivity of the collective energy to nuclear compressibility 
might be expected, as the latter generally brings less change to nuclear collisions than wide 
variations of impact parameter. In fact, e.g. in the head-on 400 MeV/nucleon Au + Au 
reaction, the net collective energy per nucleon is found the same, within ~ 1 MeV, for 
K = 200 MeV and K = 380 MeV. On the other hand, the division of kinetic energy 
into collective and excitation energy should be sensitive to particle interaction cross-sections 



which directly affect the freeze-out, see text below Eq. (21). In the aforementioned Au + Au 
reaction, the reduction of interaction cross-sections throughout the system by 30%, and the 
reduction following a parametrization |HJ ct™n = cr^w exp(— z/p/p ), with v = 0.3, lead 



to a drop in the collective energy by 16%, and 10%, respectively, compared to a standard 
calculation. 

On the other hand, a sensitivity to the nuclear compressibility can be expected for the 
anisotropy of collective energy which is associated with a delay in the start of longitudinal 
expansion compared to transverse. The delay time is of the order of 

2R-W (2R-W)( Pl - loPo ) 
7o(«o + «i) 7(A)Pi 

where W is shock width in terms of a distance in normal matter, introduced in an attempt 

to account for the mean-free-path effects, and vf 1 is shock speed with respect to compressed 

matter. Relative to the characteristic time scale for transverse expansion t exp = R/c s (char- 
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acteristic generally [f)0| , [5"l"|] , beyond the consideration in Sec. |I[B]), where c s refers to com- 
pressed matter, the delay time is 

«* = y- w i R \: : . (56) 

texp 7o(w + W ) 

Small changes in the relative delay time lead to finite changes in the energy components, 
affecting anisotropy. At any one intermediate beam energy, the ratio (|56"D is typically larger 
by ~ 17% for the stiff equation of state than for the soft, and, correspondingly, the difference 
between transverse energy per degree of freedom and longitudinal energy is larger for the stiff 
equation than for the soft. At intermediate energies the shock width in (|56"D may be expected, 



on the basis of Navier-Stokes equations [53], to be of a magnitude W ~ 2A, where A 



l/(crpo), and where o ~ 30 mb is an isotropic cross-section that yields the same transport 



properties as the differential NN cross-section |53| . This roughly conforms with W ~ 5 fm 
seen in Fig. |3|. For the stiff equation of state in the 400 MeV/nucleon Au + Au collision, the 
ratio ( |56| ) is then estimated at about 1 (c s ~ 0.6c for the state from RH equation), and for the 
soft equation at about 0.85 (c s ~ 0.45c). The reduction in the relative delay time with the 
change of compressibility leads to an increase of longitudinal energy at the cost of transverse 
energy by ~ 2 MeV/nucleon in the Au + Au simulation (see Fig. |9] for energy values), 
decreasing the anisotropy of the collective energy at freeze-out, Xcoii = — E\ oll ) / E}. oUl 

from Xcoii ^ 0.8 to Xcoii ^ 0.45. 

Following (0), a reduction in cross sections (or, more generally, interaction rates) com- 
pared to free space by ~ 20% may lead to a similar decrease in the delay time and in 
collective-energy anisotropy, as the change in the compressibility above. Indeed, in a simu- 
lation of the Au + Au reaction using a stiff equation of state, the reduction of cross section 
by 20% gives Xcoii — 0.40. A near-isotropy in collective energy could be expected, fol- 
lowing (0), for a 40-50% reduction in cross sections. In a simulation the collective energy 
components per degree of freedom become identical when cross sections are reduced by 40%. 

On the basis of ( |5^| ) and more generally, one might expect similar results from the 
dynamics when varying the size of a system by a factor, as when varying cross sections by 
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the same factor. Results would need, in fact, to be the same, up to a general rescaling 
of multiplicities, if nuclei were sharp-edged spheres. (It would not matter whether r.h.s. 
of Eqs. ([!]) and (^) were multiplied by some factor, or l.h.s. were divided by the same 
factor; see also pl| .) On the basis of the tests above, one might specifically expect that 
an isotropy of collective energy would be obtained if nuclear radius were reduced by 40%, 
i.e. mass reduced to A ~ 40. Simulation of a head-on 400 MeV/nucleon Ca + Ca reaction 
yields a result in a disagreement with this expectation, Xcoii — —0.25, demonstrating the 
importance of corona effects for a light system as compared to heavy f66|, given a similar 
width of the nuclear surface region. 

The corona effects play a lesser role when only directions pointing away from the beam 
axis are examined. Consequently, scalings expected on the basis of ( |56"D may be followed 
to a larger extent at a finite b/b max , than at b = 0. The ratio (|56|), without any modifica- 
tions, is relevant for dynamics at b/b max < 1/3. At b/b max ~ 0.25 and at the beam energy 
of 400 MeV/nucleon, the anisotropy of particle energies \ = (E y )/(E X ) -1, generally some- 
what lower than Xcoii, is X — 0-65 in a Au + Au simulation when using stiff equation of 
state, and x — 0.40, 0.45, and 0.40, respectively, when using soft equation, when reducing 
collision rates by 20%, and when turning to a Nb + Nb system with nuclei of a ~ 20% lesser 
radius than Au. 

The 400 MeV/nucleon data |67|] displayed in Fig. ^(b) rule out the soft equation of 
state with momentum dependence in the optical potential missing, in a combination with 
free-space or lesser cross-sections. After modifing the elements of calculated kinetic-energy 
tensor on the account of the fluctuations of estimated reaction-plane direction in the Au + Au 
reaction, the anisotropy of particle energies for the soft equation of state would continue to 
be about 0.6 of the anistropy for the stiff equation as above. A corresponding maximum 
value of R\ ~ 1.26 for the soft equation, as a function of the centrality, would fall below the 
data. On the other hand, the data would allow for the stiff equation of state in a combination 
with cross sections reduced by 10% at the particular beam energy compared to free space. 
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B. Hollow Structures 



In transport simulations of head-on Mo + Mo collisions at 60 and 100 MeV/nucleon, 



Moretto et al. [30] noted a formation of disk structures perpendicular to beam axis, which 
fragmented with time. Further analyses |3l,[32| revealed openings at disk centers. In the 
simulations of central collisions of heavy nuclei at 30 MeV/nucleon, bubble structures 
were observed. Hollows were attributed in 13TJ to a rarefaction wave reaching the center 
of a system, while in |33| the issues of Coulomb stability were stressed. The structures 
generated some experimental interest. Below, the mechanism of the formation of hollow 
structures is critically reassessed; the particular issues might not be of a grave importance, 
but the different arguments are being repeated in the literature. 

Despite of some increase in their width J53|, shock-like interfaces continue to develop in 
the simulations of central reactions of heavy nuclei, down to beam energies per nucleon of 
few tens of MeV. Differences between maximum density in the simulations and the RH ex- 
pectations at the lower beam energies in Fig. [|, are partially associated with a Coulomb 



repulsion [p^] . Between the shocks in reactions, an expansion develops in the directions 
parallel to shock surfaces, what gives a planarity to transient structures. At energies less 
than 35 MeV/nucleon, or for cross-sections reduced as compared to free space, and/or for 
lighter nuclei, the shock width becomes as large as the nuclear diameter and a reacting 
system equilibrates then whole. 

The onset of transient structures is marked in Fig. |6] with some increase of density in the 
vicinity of last particle collisions, towards low beam energies. The fact that structures persist 
in simulations up to the energies in excess of 100 MeV/nucleon is necessarily conditioned 
on a massive emission of rapid particles early on in reactions (if few particles were emitted, 
then structure formation would be limited to the beam energies less than 4 times the binding 
energy, i.e. ~ 30 MeV/nucleon). The simulations of a Mo + Mo system with a total mass 
of 196, within the present model, e.g. yield residues of a mass ~ 56, ~ 105, and ~ 177, 
at the beam energies of 110, 60, and 20 MeV/nucleon, respectively. Figure displays 
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baryon density, radial velocity, and entropy along and perpendicular to the beam axis from 
the 60 MeV/nucleon 6 = Mo + Mo simulations with and without Coulomb interactions. 
Emission is evidenced in extended density tails at times t > 50 fm/c. 

The argumentation in implicitly attributes hollow formation to details in a density 
distribution after equilibration, and specifically to a sharp-edged surface. Following the ar- 
guments there, one would find it difficult to understand why hollows do not develop within 
the simple solution from Sec. [Tl [B|, or within solutions in Refs. []5li|52| , or in reaction simula- 
tions at Ebeam > 150 MeV/nucleon, or, with reference to everyday experience, why rings do 
not form when spattering water on a tabletop or a wall. An edge in the density distribution 
might be built onto such a solution as in Sec. [II [B| as a perturbation. A perturbation would 



satisfy a classical wave equation. With time the initial pulse would move inwards at a speed 
of sound relative to remaining matter, and emerge from the center as a wiggle in the case of 
a cyllindrical symmetry, with a peak in front followed by a dip, and as as a dip in the case 
of a spherical symmetry. (Phase shift at the center for the Fourier components of the pulse 
is 7r/2 in the case of a cyllindrical symmetry and 7r in the case of a spherical symmetry.) 
As can be seen in Fig. the expansion in simulations develops past the maximum in dv/dr 
in transverse direction, at t ~ 36 fm/c, without any noticable depletion in the central den- 
sity relative to outer regions. This is due to a continuous behavior of parameters of nuclear 
matter after equilibration. A depletion appears later and is associated with the changes in 
outer regions of an expanding system. 

As the nuclear system expands at a low entropy S/A ~ 1, see Fig. [32], down to subnormal 
densities, a negative pressure develops, cf. Fig. |33|. As a result, outer edges of the expanding 
matter begin to deccelerate (the vacuum pressure is zero); see the plateau in transverse 
velocity in Fig. ^ developing after central density falls below 0.9 po- in consequence, with 
time the density becomes enhanced at the outer edges as compared to the center. As the 
system is driven in the meantime into a region of adiabatic instability within thermodynamic 
parameters (see also | j51)| , |9"2"| |), characterized by {dP/dp)s/A < 0, the pressure becomes more 
negative at the outskirts with a higher density, than at the center. From the center then 
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the matter continues to move out and accumulate in the outer regions. The system is being 
cooled off all the time by emission. Until t ~ 70 fm/c from the Mo + Mo system displayed 
in Fig. |32], particles with a total mass of ~ 35 are emitted. At time t ~ 100 fm/c the 
total mass of emitted particles rises to ~ 60. Apart from the pressure within the system, 
effects of which are local, Coulomb interactions act. With time, these render a stability to 
a hollow structure that forms as the matter contracts at the outskirts of the system. When 
the Coulomb interactions are missing or are weak, the matter having first contracted at the 



outskirts, gradually collapses onto the center, see Fig. |32](b). This is the case in simulations 
of such light systems as Ca + Ca. 

Deficiency of the single-particle models is the lack of fluctuations. Because of the fluctu- 
ations, a system might, in reality, not condense into a ring or a bubble, but into a number 
of clusters. If a structure were formed, it would later break, due to fluctuations, into clus- 
ters of which diameter would be likely commensurate with the minor diameter of a ring, 
or the thickness of a bubble. In any of these cases a large number of IMF would be pro- 
duced, slow in the c.m.s. on a scale of collective velocities characterizing the early emission. 
The total mass of these IMF, e.g. at the beam energy above, could near half of the mass 
of the system. The unusual abundance of these fragments, already indicated in 153], might 



thus serve as a signature of the development of an adiabatic instability over a macroscopic 
region in reactions. The latter underlies both of the scenarios outlined above. Evidence for 
a planarity in the emission of the slow IMF would additionally show that the instability sets 
in on a dynamic path for the system. Search for signatures of transient structures in light 
particles might, on the other hand, be hopeless. These particles are primarily emitted when 
the fate of matter left behind has not yet been decided. 

VIII. CONCLUSIONS 

Dynamics of energetic symmetric heavy-ion reactions has been analyzed at a qualitative 
level, with an attention directed at the collective behavior of nuclear matter. The analysis 
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relied on transport-model reaction-simulations, analytic and near-analytic considerations, 
and measurements. 

In the simulations of head-on reactions, nearly completely developed shock-fronts may 
be identified, propagating into the projectile and target, and separating the hot matter 
from normal. In-between the shocks, at finite impact parameters, a tangential discontinuity 
develops, that spreads with time. Because of small nuclear sizes and finite widths of the 
discontinuities, they strongly overlap. Hot matter exposed to the vacuum in sidewards 
directions begins to expand into these directions. Expansion in the shock direction is delayed 
and does not acquire same strength. The collective expansion affects angular distributions, 
mean-energy components, shapes of spectra and mean energies of different particles emitted 
into any one direction, and particle yields. 

The stronger expansion perpendicular to the shock direction leads to a cm. 90° peaking 
in the polar angle in head-on collisions of heavy nuclei, and to a peaking in the azimuthal 
angle out of the reaction plane in semicentral collisions. These effects contrast with the 
naive expectations regarding shocks, where one expects an enhancement of the emission in 
shock direction. Unambigous experimental evidence exists only for peaking of the emission 
in azimuthal angle. Evidence for a motion of matter behind shocks, as characteristic for 
the weak discontinuity, exists in the strength of the sideward flow in collisions. The strong 
overlap of the discontinuities dampens this motion, nonetheless, to the extent that spectator 
pieces can emerge in the high-energy collisions. 

The maximum density in reactions falls somewhat below the RH expectations. As a func- 
tion of the impact parameter, this density follows well the scaling relation ( p6|) , falling off 
quadratically with the impact parameter, while the entropy at the center gets close to that 
characteristic for stopped matter. It follows then that the thermodynamic parameters of 
matter in the central region of reactions can be regulated by changing the beam energy 
and selecting impact parameter. With regard to expansion, its relative strength along and 
perpendicular to the shock front motion depends on the magnitude of time necessary for 
a shock propagation through a nucleus, relative to the time necessary for the development 
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of expansion. The relative delay time is larger for the stiff equation of state than for the 
soft, and the anisotropies in emission, as quantified e.g. using collective-energy or particle- 
energy components, are larger for the stiff equation of state than for the soft. The delay 
time depends, moreover, on the interaction cross sections that determine the width of the 
shock front and the size of a region within which original equilibration occurs. The division 
of the energy into collective and thermal for a given system at low b, occuring at freeze-out, 
depends solely on the cross sections. 

The effects of expansion on particle spectra at any one cm. angle in central reactions, 
are the slopes decreasing with increasing particle mass and downward concavities in baryon 
spectra. The effect on pion spectra is stronger than following from naive nonrelativistic 
considerations. A pion spectrum is overall flatter than could be expected from latter con- 
siderations and an upward concavity can develop. Resonance decays may, however, fill up 
this concavity. With regard to pion yields, the missing energy considered by Stock et al. 
when analysing data, may be identified with the collective energy not available locally for 
pion production. The number pions and deltas actually decreases within a reacting system 
towards freeze-out, as the collective motion with participant region acquires strength. Over- 
all, the effect of expansion reduces pion yields by ~ 50% compared to the situation without 
expansion. 

Collective expansion is not fully isentropic towards freeze-out, with produced entropy 
being associated with the heat conduction. Besides, in reactions the entropy is produced in 
shock fronts and within a weak discontinuity on account of the shear viscosity. The shock 
contribution to the net entropy maximizes in head-on reactions. In low-energy collisions, 
collective expansion initiates the formation of hollow structures. The process further involves 
negative pressure and adiabatic instability, and finally Coulomb interactions. 

As may be apparent, a wealth of physical phenomena is associated with the collective 
behavior of matter in central collisions. This is in stark contrast with early conclusions on 
heavy-ion collisions based on limited inclusive data which, dominated by peripheral colli- 
sions, appeared consistent with a featureless global equilibrium within the participant region. 
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FIGURES 

FIG. 1. Proton (solid lines) and neutron (dashed lines) density profiles from solving the TF 
equations, together with the empirical charge density profiles (dotted lines) |4j| for 40 Ca and 208 Pb. 

FIG. 2. Contour plots of baryon density in the reaction plane in Au + Au collisions at 400 
MeV/nucleon. The displayed contour lines are for the densities p/po = 0.1, 0.5, 1, 1.5, and 2. 

FIG. 3. Baryon density (top panels), radial velocity (center panels), and entropy per baryon 
along (solid lines) and perpendicular (dashed lines) to the beam axis at different indicated times, 
in the b = collision at a beam energy of 400 MeV/nucleon. 

FIG. 4. Lines show the baryon density (top panel), temperature (center panel), and total 
entropy per baryon (bottom panel) expected behind a developed shock front at b = 0, as a function 
of beam energy, from Eq. (|l6|). Filled circles indicate the maximum density from the simulations of 
6 = Au + Au collisions at different beam energies. Crosses indicate the maximum density, scaled 
by a 7- factor according to Eq. (^), from the simulations of 1 GeV/nucleon Au + Au collisions at b 
= 3, 6, 8, 9, 10, 11, and 12 fm, plotted against the equivalent beam energy following from Eq. (|26|). 

FIG. 5. Time dependence of the baryon density at r = (top panel) and the participant 
transverse collective energy per nucleon (bottom panel) in the b = Au + Au reaction at 400 
MeV/nucleon. 

FIG. 6. Mean baryon densities in units of normal density in the vicinity of last collisions of 
emitted particles from the head-on Au + Au reactions (solid line), and in the vicinity of collisions 
from which the final light-clusters originate (dashed line) , as a function of beam energy. 

FIG. 7. Distribution in the spherical angle of protons (circles), deuterons (triangles), tritons 
(inverted triangles), pions (diamonds), and pions emitted within first 14 fm/c (stars), as a function 
of polar angle from the simulation of a b = Au + Au collision at 1 GeV/nucleon. The solid line 
indicates the distribution of tritons expected on the basis of the formula (|23|). 
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FIG. 8. Transverse and longitudinal components of the mean energies of particles emitted 
from 6 = Au + Au reactions, as a function of beam energy. 

FIG. 9. Transverse and longitudinal components of the mean collective energy per nucleon 
at the time of last collision of baryons stemming from 6 = Au + Au reactions, as a function of 
beam energy. 

FIG. 10. The initial discontinuity between the projectile and target velocities breaks at a 
finite 6 into two shock fronts propagating into the projectile P and target T (thick solid lines), and 
a weak tangential discontinuity in-between (dotted line). Shock- front position within projectile at 
a later time is indicated with a dashed line. 

FIG. 11. Baryon density p (top panel), normal and tangential velocity-components v n and 
vt to the discontinuities in nuclear- matter (center panel), coinciding with radial and polar veloc- 
ity-components, and entropy per baryon (bottom panel) as a function of the distance from the 
center of a 400 MeV/nucleon Au + Au system at 6 = 7 fm, along the normal to the discontinuities. 
At the given time t = 13.5 fm/c the discontinuities (at their centers) are inclined at an angle 
a ~ 52° relative to the beam axis. 

FIG. 12. Mean energy component out of the reaction plane at y = (top panel), and at 
all rapidities (center panel), and a lower of the eigenvalues of a relativistic tensor (s*- 7 ), associated 
with a direction within the reaction plane (bottom panel), for protons (circles), deuterons (squares), 
tritons (triangles), and helions (diamonds), emitted from 400 MeV/nucleon Au + Au reactions, as 
a function of the impact parameter. 

FIG. 13. Mean transverse momentum component within the reaction plane, divided by mass, 
as a function of rapidity for protons (circles), deuterons (triangles), tritons (inverted triangles), 
and pions (diamonds) from 1 GeV/nucleon Au + Au collisions at 6 = 3 fm (top panel) and 6 = 7 
fm (bottom) panel. For collisions at 6 = 7 fm, also the mean momentum for pions emitted within 
first 16 fm/c is shown (stars). 
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FIG. 14. Left panel shows the flow parameter ( |28|) in Au + Au reactions at 250 and 400 
MeV/nucleons, from the measurements of Ref. |3l|], as a function of the reduced participant proton 
multiplicity. Right panel shows the flow parameter from the simulations of Au + Au reactions at 
400 MeV/nucleon, as a function of impact parameter. 

FIG. 15. Parameters associated with the kinetic-energy tensor in symmetric reactions at 400 
MeV/nucleon. (a) Flow angle from the calculated mean kinetic-energy tensor, as a function of 
reduced impact parameter in Au + Au (stars) and Ca + Ca (diamonds) reactions, (b) Ratio of the 
out-of-reaction-plane matrix-element to the lower in-plane eigenvalue (crosses), as a function of the 



normalized participant proton multiplicity in Au + Au reactions, from measurements of Ref. [67], 
compared to the ratio of eigenvalues (stars) from a calculated tensor with matrix elements modified 
on account of fluctuations of the estimated reaction plane about the true plane. 

FIG. 16. Mean kinetic energy divided by temperature in thermal equilibrium, as a function 
of temperature divided by mass. 

FIG. 17. Baryon density along (solid lines) and perpendicular (dashed lines) to the beam axis 
at different indicated times, in the b = La + La collision at a beam energy of 800 MeV/nucleon. 

FIG. 18. Time dependence of baryon density at r = (top panel), participant collective energy 
per nucleon (center panel), and total entropy in the b = La + La reaction at 800 MeV/nucleon 
from the standard calculation (solid line) and from the calculation with an interchange of particle 
positions (short-dashed lines). 

FIG. 19. Time dependence of elastic NN collision-rate (top panel), rates of A production in 
NN collisions and of absorption in AN collisions (center panel), and rates of A formation in 7rN 
collisions and of A decay (bottom panel). 
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FIG. 20. Time dependence of the number of pions and deltas normalized using the baryon 
number, in the 6 = La+La reaction at 800 MeV/nucleon, from a standard calculation (solid 
line) and from a calculation with an interchange of particle positions (short-dashed line). The two 
horizontal long-dashed lines show the equilibrium number of pions and deltas for a baryon density 
p = 2.25/00 an d an energy per baryon of 1133 MeV and 1103 MeV, respectively. Dots show the 
number obtained by integrating over space the local equilibrium values of pion and delta density 
in a standard calculation. 

FIG. 21. Isentropes in the plane of the ratio of pion-and-delta density to baryon density vs 
the baryon density (top panel), and in the plane of the temperature vs the baryon density (bottom 
panel). Approximate locations of the 800 MeV/nucleon 6 = La+La system at t = 14fm/c in the 
two planes are marked with crosses. Dashed lines show a portion of a chemical equilibrium path 
with regard to the pion and delta number. 

FIG. 22. Momentum distribution of protons (circles), deuterons (triangles), helions (inverted 
triangles), negative (squares) and neutral (diamonds) pions from central 800 MeV/nucleon La + La 
reaction, in the vicinity of 90° in the cm. Left panel shows the results of calculations at 6 = 1 fm 
(filled symbols) and data of Ref. J7(| (open symbols). Right panel shows the results of calculations 
with particle positions interchanged during evolution. Solid lines in the left panel indicate results 
of the instantaneous freeze-out model, for baryons and neutral pions. Long- and short-dashed 
lines for neutral pions indicate a contribution of free pions at freeze-out and a contribution from 
A decays, respectively. Straight parallel lines in the right panel serve to guide the eye. 

FIG. 23. Ratio of the mean pion multiplicity to the number of participant nucleons in central 
symmetric reactions, as a function of laboratory energy. The solid and dashed lines represent, 
respectively, the results of standard calculations and of calculations with particle positions inter- 
changed, for the La+La 6 = reactions. The circles represent data of Ref. ||. 
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FIG. 24. Distribution of collective velocities for baryons from 800 MeV/nucleon La + La reac- 
tion at b = 1 fm, calculated in the vicinity of last collisions, within an angular range of velocities 
9 = 60 - 120°. 

FIG. 25. Negative-pion distribution in the vicinity of cm. 90° in the central 1.8 GeV/nucleon 
Ar + KC1 reaction. Insets show pion asymmetry-coefficient defined in Eq. (^). Open circles in 
both panels indicate data of Ref. ||; distributions are normalized using event cross-section given 
there. Filled circles in the left panel indicate results of the present model at b = 1.7 fm. Squares 
and diamonds in the right panel indicate, respectively, the results obtained using the cascade code 



[73] with anisotropic and isotropic N + N «-> N + A cross-sections. It may be mentioned that the 



present model yields results indistinguishable from those from the cascade code [73], when optical 
potential is switched off and same assumptions on scattering are adopted. 

FIG. 26. Distribution of neutral pions with cm. momenta p > 200 MeV/c (top panel) and 
protons (bottom panel) in the azimuthal angle <p' with respect to the reaction plane, at a momentum 
component within the reaction plane p z ' = 0. The angle G is the angle of rotation of the p z '-axis 
from the p z -axis in the direction of p x . 

FIG. 27. Contour lines (solid) for baryon density in the reaction plane in 1 GeV/nucleon Au 
+ Au collision at b = 7 fm and t = 15 fm/c, with overlaid contour lines (dashed) for pion and delta 
density. Baryon contour lines are for density values p/po = 0.1, 0.5, 1, 1.5, and 2. Pion and delta 
lines are for density values (p n + pa)/ Po = 0.02, 0.1, and 0.2. 

FIG. 28. Left panel shows entropy per nucleon as a function of reduced participant proton 
multiplicity, determined from 400 MeV/nucleon Au + Au wide-angle data of Ref. [p0[| . Right panel 
shows, as a function of impact parameter, the entropy per nucleon associated with nucleons emitted 
into wide angles from reaction simulations (circles), cf. |f35]| , and predictions for the entropy from 
Eq. ( |53"| ) (dashed line) and from Eq. ( |54] ) (solid line). 
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FIG. 29. Ratios of cluster yields to proton yields at wide angles, cf. [35], as a function of bom- 
barding energy in the head-on Au+Au reactions. The solid and dashed lines represent, respectively, 
the results from standard calculations and from calculations with particle positions interchanged. 
The dots represent data of Ref. |2C|. 



FIG. 30. Distribution of entropy per baryon in the vicinity of last collision points in the 400 
MeV/nucleon 6 = Au + Au reaction simulation. 

FIG. 31. Impact-parameter dependence of the mean participant collective energy per nucleon 
(stars) and transverse collective energy per nucleon (diamonds), calculated as in Eq. (|2l|) , in 
Au + Au reactions at 400 MeV/nucleon. 

FIG. 32. Baryon density (top panels), radial velocity (center panels), and entropy per baryon 
along (solid lines) and perpendicular (dashed lines) to the beam axis at different indicated times, 
in the head-on Mo + Mo reaction at 60 MeV/nucleon, from a simulation with (a) and without (b) 
Coulomb interactions. 

FIG. 33. Isentropes in the plane of pressure vs density for the stiff equation of state. 
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